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Resume 


Dans  ce  Rapport,  les  grandes  lignes  d'un  modele  analytique  permettant  le  calcul  des  niveaux  de  gris  delivres  par  une 
camera  CCD  lors  des  prises  de  vues  par  la  methode  de  Velodmetrie  par  Images  de  Particules  (PIV)  sont  rappelees. 
L'image  des  particules  eclairees  par  le  plan  laser  est  determinee  a  partir  d'une  relation  de  convolution  entre  la  reponse 
impulsionnelle  de  I'objectif  et  l'image  geometrique  de  ces  particules,  et  la  camera  CCD  est  modelisee  par  ses  effets 
d'integration  spatiale,  de  gain  et  de  quantification.  Ce  modele  est  valide  sur  un  cas  precis  de  prises  de  vues  reelles 
de  PIV. 

Les  principes  des  methodes  classiques  de  correlation  (autocorrelation  et  correlation  croisee)  et  d'interpolation  sub¬ 
pixel  (barycentre,  ajustement  par  fonction  parabolique  et  gaussienne)  permettant  la  determination  du  deplacement 
des  particules  sont  egalement  rappeles. 

Ce  Rapport  se  poursuit  par  des  simulations  numeriques,  basees  sur  le  modele  precedent,  correspondant  a  deux 
configurations  de  particules  :  la  premiere  configuration  concerne  le  cas  ideal  d'une  seule  particule  par  aire 
d'interrogation  et  la  seconde  configuration,  le  cas  plus  realiste,  d'un  nombre  variable  de  particules  par  aire 
d'interrogation.  Dans  ces  deux  configurations,  les  sources  de  bruit  ne  sont  pas  prises  en  compte.  Ces  simulations 
fournissent  des  resultats  statistiques  importants  concernant  les  erreurs  commises  par  |es  methodes  classiques  de 
correlation  et  d'interpolation  sub-pixel  lors  de  la  determination  de  deplacement  des  particules.  De  plus,  ces 
simulations  mettent  en  evidence  les  faiblesses  des  methodes  employees  (production  de  "faux  vecteurs").  Enfin,  ce 
Rapport  aborde  le  cas  de  configurations  de  particules  similaires  a  celles  rencontrees  sur  des  prises  de  vues  reelles  de 
PIV  (nombre  variable  de  particules,  bruit,  gradient  de  vitesses). 


Mots-des 

Autocorrelation 
Camera  CCD 
Correlation  croisee 
Erreurs  de  deplacement 
Faux  vecteurs 
Images  de  particules 


Interpolation  sub-pixel 

Mecanique  des  fluides 

Methodes  de  mesure 

Optique  de  Fourier 

Simulations  numeriques 

Velodmetrie  par  images  de  particules  (PIV) 
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Zusammenfassung 


In  diesem  Bericht  wird  ein  analytisches  Modell  zur  Berechnung  der  Grauwerte  einer  CCD-Kamera  bei  Aufnahmen  mit- 
tels  der  PIV-Methode  (Particle  Image  Velocimetry)  umrissen.  Das  Bild  der  durch  den  Lichtschnitt  beleuchteten  Partikel 
wird  ausgehend  von  einer  Faltungsbeziehung  zwischen  der  Impulsantwort  des  Objektivs  und  dem  geometrischen  Bild 
dieser  Partikel  bestimmt,  und  die  CCD-Kamera  wird  anhand  ihrer  Wirkungen  betreffend  raumliche  Integration,  Ver- 
starkung  und  Quantifizierung  modelliert.  Dieses  Modell  wird  an  einem  prdzisen  Fall  realer  PIV-Aufnahmen  validiert. 
Die  Prinzipien  der  klassischen  Methoden  der  Korrelation  (Autokorrelation  und  Kreuzkorrelation)  und  der  Subpixel- 
Interpolation  (Schwerpunkt,  Anpassung  durch  parabolische  und  GauBsche  Funktion),  die  die  Bestimmung  der  Parti- 
kelverschiebung  ermoglichen,  werden  ebenfalls  nochmals  erlautert. 

Des  Weiteren  enthdlt  dieser  Bericht  numerische  Simulationen  auf  der  Basis  des  vorstehend  angefuhrten  Modells,  ent- 
sprechend  zwei  Partikelkonfigurationen:  die  erste  Konfiguration  betrifft  den  idealen  Fall  eines  einzigen  Partikels  pro 
Abfragebereich  und  die  zweite  Konfiguration  den  realistischeren  Fall  einer  variablen  Partikelzahl  pro  Abfragebereich. 
Bei  beiden  Konfigurationen  werden  die  Rauschquellen  nicht  berucksichtigt. 

Diese  Simulationen  liefern  wichtige  statistische  Ergebnisse  uber  die  Fehler  durch  die  klassischen  Methoden  der  Korre¬ 
lation  und  Subpixel-Interpolation  bei  der  Bestimmung  der  Partikelverschiebung.  AuBerdem  zeigen  diese  Simulationen 
die  Schwachstellen  der  verwendeten  Methoden  auf  (Produktion  von  "falschen  Vektoren"). 

SchlieBlich  wird  in  diesem  Bericht  der  Fall  von  Partikelkonfigurationen,  ahnlich  denjenigen  bei  realen  PIV-Aufnahmen, 
angesprochen  (variable  Partikelzahl,  Rauschen,  Geschwindigkeitsgradient). 


Schlagworter 
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Verschiebungsfehler 
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Zeitpunkt  t]  und  der  Kreis  zum  Zeitpunkt  t2 . 44 

Fig.  1  9  -  Histogramme  de  repartition  du  nombre  de  couples  de  particules  presents 
dans  la  fenetre  d'interrogation  aux  instants  t]  et  t2 

Saulendiagramm  der  Verfeilung  derZahl  von  Parfikelpaaren ,  die  sich  zu  den 


Zeitpunkten  t]  und  f2  im  Abfragefenster  befinden . 44 

Fig.  20  -  Histogramme  de  repartition  du  deplacement  ^/ax2  +  Ay2 

Saulendiagramm  der  Verfeilung  der  Verschiebung  ^Ax2  +  Ay2  . 45 


Fig.  21  a,b  -  Convergence  de  Perreur  de  deplacement  (en  pixel)  obtenue  par  la  methode 

du  centre  des  masses  en  autocorrelation  (fig.  21a)  et  en  correlation  croisee  (fig.  21b) 
Konvergenz  des  Verschiebungsfehlers  (in  Pixel)  aufgrund  der  Methode  des  Massen- 
mittelpunkts  bei  der  Autokorrelation  (Fig.  21  a)  und  der  Kreuzkorrelation  (Fig.  21b) . 45 

Fig.  21  c,d  -  Convergence  de  Perreur  de  deplacement  (en  pixel)  obtenue  par  la  methode  de 

I'ajustement  parabolique  en  autocorrelation  (fig.  21c)  et  en  correlation  croisee  (fig.  21  d) 
Konvergenz  des  Verschiebungsfehlers  (in  Pixel)  aufgrund  der  Methode  der  paraboli- 
schen  Anpassung  bei  der  Autokorrelation  (Fig.  21  c)  und  der  Kreuzkorrelation  (Fig.  2  Id)..  45 
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Fig.  21  e,f  -  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenue  par  la  methode 

de  Fajustement  gaussien  en  autocorrelation  (fig.  21e)  et  en  correlation  croisee  (fig.  2 If) 
Konvergenz  des  Verschiebungsfehlers  (in  Pixel)  aufgrund  der  Methode  der  GauBschen 
Anpassung  bei  der  Autokorrelation  (Fig.  21  e)  und  der  Kreuzkorrelation  (Fig.  2 If)  . 46 

Fig.  22  a,b  -  Repartition  des  erreurs  de  deplacement  (en  pixel)  obtenues  par  la  methode  du 
centre  des  masses  en  autocorrelation  (fig.  22a)  et  en  correlation  croisee  (fig.  22b) 

Verteilung  der  Verschiebungsfehler  (in  Pixel)  aufgrund  der  Methode  des  Massen- 
mittelpunkts  bei  der  Autokorrelation  (Fig.  22a)  und  der  Kreuzkorrelation  (Fig.  22b) . 46 

Fig.  22  c,d  -  Repartition  des  erreurs  de  deplacement  (en  pixel)  obtenues  par  la  methode  de 

Fajustement  parabolique  en  autocorrelation  (fig.  22c)  et  en  correlation  croisee  (fig.  22d). 
Verteilung  der  Verschiebungsfehler  (in  Pixel)  aufgrund  der  Methode  der  parabolischen 
Anpassung  bei  der  Autokorrelation  (Fig.  22c)  und  der  Kreuzkorrelation  (Fig.  22d)  . 47 

Fig.  22  e,f  -  Repartition  des  erreurs  de  deplacement  (en  pixel)  obtenues  par  la  methode 

de  lajustement  gaussien  en  autocorrelation  (fig.  22e)  et  en  correlation  croisee  (fig.  22f) 
Verteilung  der  Verschiebungsfehler  (in  Pixel)  aufgrund  der  Methode  der  GauBschen 


Anpassung  bei  der  Autokorrelation  (Fig.  22e)  und  der  Kreuzkorrelation  (Fig.  22f)  . 47 

Fig.  23  -  Repartition  des  particules  donnant  lieu  a  un  "faux  vecteur"  :  instant  t-j  (croix) 
et  instant  t2  (rond) 

Verteilung  der  zu  einem  "falschen  Vektor "  fuhrenden  Partikel :  Zeitpunkt  t j 

(Sternchen)  und  Zeitpunkt  t2  (Kreis) . 48 

Fig.  24  -  Image  restituee  par  la  camera  CCD  (aire  de  32  x  32  pixels) 

Durch  die  CCD-Kamera  wiedergegebenes  Bild  (Bereich  von  32  x  32  Pixel) . 48 

Fig.  25  -  Repartition  tridimensionnelle  des  niveaux  de  gris  correspondant  aux 
figures  23  et  24 

Dreidimensionale  Verteilung  der  Grauwerte  entsprechend  den  Figuren  23  und  24 . 48 


Fig.  26  -  Ajustement,  par  des  fonctions  parabolique  et  gaussienne,  des  histogrammes 
du  pic  principal  de  la  figure  d'autocorrelation  correspondant  a  la  figure  25  ; 
a)  histogramme  suivant  Faxe  x  ;  b)  histogramme  suivant  Faxe  y 
Anpassung  -  durch  parabolische  und  GauBsche  Funktion  -  der  Saulendiagramme 
der  Ftauptspitze  der  Autokorrelationsabbildung  entsprechend  der  Figur  25; 
a)  Saulendiagramm  gemaB  der  x-Achse;  b)  Saulendiagramm  gemaB  der  y-Achse . . 48 

Fig.  27  -  Ajustement,  par  des  fonctions  parabolique  et  gaussienne,  des  histogrammes 

du  pic  principal  de  la  figure  de  correlation  croisee  correspondant  a  la  figure  25  ; 

a)  histogramme  suivant  Faxe  x  ;  b)  histogramme  suivant  Faxe  y 

Anpassung  -  durch  parabolische  und  GauBsche  Funktion  -  der  Saulendiagramme 

der  Flauptspitze  der  Kreuzkorrelationsabbildung  entsprechend  der  Figur  25; 

a)  Saulendiagramm  gemaB  der  x-Achse;  b)  Saulendiagramm  gemaB  der  y-Achse . 49 

Fig.  28  a,b  -  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenues  apres 
elimination  des  "faux  vecteurs"  par  la  methode  du  centre  des  masses  en 
autocorrelation  (fig.  28a)  et  en  correlation  croisee  (fig.  28b) 

Konvergenz  des  Verschiebungsfehlers  (in  Pixel)  nach  Beseitigung  der  "falschen 

Vektoren 11  durch  die  Methode  des  Massenmittelpunkts  bei  der  Autokorrelation 

(Fig.  28  a)  und  der  Kreuzkorrelation  (Fig.  28b) . 50 

Fig.  28  c,d  -  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenues  apres 

elimination  des  "faux  vecteurs"  par  la  methode  d'ajustement  parfonction  parabolique 
en  autocorrelation  (fig.  28c)  et  en  correlation  croisee  (fig.  28d) 

Konvergenz  des  Verschiebungsfehlers  (in  Pixel)  nach  Beseitigung  der  "falschen  Vektoren" 

durch  die  Methode  der  Anpassung  durch  parabolische  Funktion  bei  der 

Autokorrelation  (Fig.  28c)  und  der  Kreuzkorrelation  (Fig.  28d)  . 51 
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Fig.  28  e,f  -  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenues  apres  elimination 
des  "faux  vecteurs"  par  la  methode  d'a|ustement  par  fonction  gaussienne 
en  autocorrelation  (fig.  28e)  et  en  correlation  croisee  (fig.  28f) 

Konvergenz  des  Verschiebungsfehlers  (in  Pixel)  nach  Beseitigung  der  11 falschen  Vekforen " 
durch  die  Methode  der  Anpassung  durch  GauBsche  Funkfion  bei  der  Autokorrelafion 


(Fig.  28e)  und  der  Kreuzkorrelafion  (Fig.  28f)  . 51 

Fig.  29  a,b  -  Pourcentages  de  simulations  dont  I'erreur  est  superieure  a  E0  ; 
a)  autocorrelation  ;  b)  correlation  croisee 
Prozentsatz  von  Simulafionen  mil  einem  Fehler  groBer  als  EQ ; 

a)  Autokorrelafion;  b)  Kreuzkorrelafion . 52 

Fig.  29  c,d  -  Pourcentages  de  simulations  dont  I'erreur  est  superieure  a  E0  ; 
c)  autocorrelation  ;  d)  correlation  croisee 
Prozentsatz  von  Simulafionen  mit  einem  Fehler  groBer  als  E0 ; 

c)  Autokorrelafion ;  d)  Kreuzkorrelafion . 52 

Fig.  30a  -  Erreur  de  deplacement  en  fonction  du  deplacement  fractionnaire 


(cas  de  I'autocorrelation)  :  moyenne  (en  pixel)  de  I'erreur  de  deplacement 

Verschiebungsfehler  in  Abhangigkeit  von  der  anfeiligen  Verschiebung 

(Fall  der  Autokorrelafion):  mifflerer  Verschiebungsfehler  (in  Pixel) . 53 

Fig.  30b  -  Erreur  de  deplacement  en  fonction  du  deplacement  fractionnaire 

(cas  de  I'autocorrelation)  :  ecart  type  (en  pixel)  de  I'erreur  de  deplacement 
Verschiebungsfehler  in  Abhangigkeit  von  der  anfeiligen  Verschiebung 

(Fall  der  Autokorrelafion):  Standardabweichung  (in  Pixel)  des  Verschiebungsfehlers . 53 

Fig.  30c  -  Erreur  de  deplacement  en  fonction  du  deplacement  fractionnaire 

(cas  de  la  correlation  croisee)  :  moyenne  (en  pixel)  de  I'erreur  de  deplacement 
Verschiebungsfehler  in  Abhangigkeit  von  der  anfeiligen  Verschiebung 

(Fall  der  Kreuzkorrelafion):  mifflerer  Verschiebungsfehler  (in  Pixel) . 53 

Fig.  30d  -  Erreur  de  deplacement  en  fonction  du  deplacement  fractionnaire 

(cas  de  la  correlation  croisee)  :  ecart  type  (en  pixel)  de  I'erreur  de  deplacement 
Verschiebungsfehler  in  Abhangigkeit  von  der  anfeiligen  Verschiebung 

(Fall  der  Kreuzkorrelafion):  Standardabweichung  (in  Pixel)  des  Verschiebungsfehlers . 53 

Fig.  31a  -  Erreur  de  deplacement  (moyenne,  en  pixel)  en  fonction  du  rayon  des  particules 
(cas  de  I'autocorrelation) 

Verschiebungsfehler  ( mifflerer ,  in  Pixel)  in  Abhangigkeit  vom  Parfikelradius 

(Fall  der  Autokorrelafion) . 54 

Fig.  31b-  Erreur  de  deplacement  (ecart  type,  en  pixel)  en  fonction  du  rayon  des  pariicules 
(cas  de  ('autocorrelation) 

Verschiebungsfehler  ( Standardabweichung ,  in  Pixel)  in  Abhangigkeit  vom  Parfikelradius 
(Fall  der  Autokorrelafion) . 54 

Fig.  31c  -  Erreur  de  deplacement  (moyenne,  en  pixel)  en  fonction  du  rayon  des  particules 
(cas  de  la  correlation  croisee) 

Verschiebungsfehler  ( mifflerer ,  in  Pixel)  in  Abhangigkeit  vom  Parfikelradius 

(Fall  der  Kreuzkorrelafion) . 54 

Fig.  31  d  -  Erreur  de  deplacement  (ecart  type,  en  pixel)  en  fonction  du  rayon  des  particules 
(cas  de  la  correlation  croisee) 

Verschiebungsfehler  ( Standardabweichung ,  in  Pixel)  in  Abhangigkeit  vom 

Parfikelradius  (Fall  der  Kreuzkorrelafion) . 54 
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Fig.  32a  -  Position  des  particules  de  rayon  egal  a  3  pixels  dans  une  aire  d'interrogation 
de  32  X  32  pixels  :  la  croix  correspond  a  la  position  des  particules  a  Pinstant  t-| 
et  le  rond  a  Pinstant  t2 

Position  der  Parfikel  mil  einem  Radius  gleich  3  Pixel  in  einem  Abfragebereich 
von  32  x  32  Pixel:  das  Sfernchen  entspricht  der  Position  der  Parfikel  zum 


Zeitpunkt  t]  und  der  Kreis  zum  Zeitpunkt  t2 . . 55 

Fig.  32b  -  Image  restituee  par  la  camera  CCD 

Durch  die  CCD-Kamera  wiedergegebenes  Bild . 55 

Fig.  33a  -  Positions  de  particules  de  rayon  0,4  pixel 

Positionen  von  Partikeln  mit  einem  Radius  von  0,4  Pixel . 55 


Fig.  33b  -  Niveaux  de  gris  correspondent  a  la  figure  33a  (les  pics  correspondent  pour 
la  piupart  a  des  superpositions  de  particules) 

Grauwerte  entsprechend  der  Figur  33a  (die  Spitzen  entsprechen  zumeist 
Partikeluberlagerungen) . 55 

Fig.  34  -  Representations  bi-  et  tridimensionnelle  des  niveaux  de  gris  d'une  aire 
d'interrogation  obtenue  lors  d'une  prise  de  vue  reelle  de  PIV  [14,15] 

Zwei-  und  dreidimensionale  Darstellung  der  Grauwerte  eines  Abfragebereichs 

bei  einer  realen  PIV-Aufnahme  [14,  15] . 56 

Fig.  35a  -  Position  des  particules  dans  Paire  d'interrogation:  la  croix  correspond 
a  la  position  des  particules  a  Pinstant  t]  et  le  rond  a  Pinstant  t2 
Position  der  Parfikel  im  Abfragebereich:  das  Sternchen  entspricht  der 


Partikelposition  zum  Zeitpunkt  tj  und  der  Kreis  zum  Zeitpunkt  ^ . 57 

Fig.  35b  -  Image  restituee  par  la  camera  CCD 

Durch  die  CCD-Kamera  wiedergegebenes  Bild . . . . 57 

Fig.  35c  -  Repartition  tridimensionnelle  des  niveaux  de  gris  correspondant  a  la  figure  35a 

Dreidimensionale  Verteilung  der  Grauwerte  entsprechend  der  Figur  35a . 57 

Fig.  36a  -  Fonction  d'autocorrelation 

Autokorrelationsfunktion . 57 

Fig.  36b  -  Fonction  de  correlation  croisee 

Kreuzkorrelationsfunktion . 57 

Fig.  37  -  Fonction  d'autocorrelation  (p(t)  d'une  fonction  constante  f(x)  sur  un  intervalle 

Autokorrelationsfunktion  cpf  t)  einer  konstanten  Funktion  f(x)  an  einem  Interval  1 . 58 


Fig.  38  -  Pourcentage  de  simulations  dont  I'erreur  de  deplacement  est  1  0  fois  inferieure 
apres  elimination  du  bruit  a  I'erreur  de  deplacement  obtenue  avant  elimination 
du  bruit 

Prozentsatz  von  Simulationen,  deren  Verschiebungsfehler  nach  der  Beseifigung 
des  Rauschens  1  Omal  kleiner  ist  als  der  Verschiebungsfehler  vor  der  Beseitigung 


des  Rauschens . 60 

Fig.  39a  -  Exemple  de  niveaux  de  gris  avant  elimination  du  bruit 

Beispiel  von  Grauwerten  vor  der  Beseitigung  des  Rauschens . 60 

Fig.  39b  -  Niveaux  de  gris  apres  elimination  du  bruit  avec  un  nombre  seuil  egal  a  0,1 

Grauwerte  nach  der  Beseitigung  des  Rauschens  mit  einem  Schwellenwert  gleich  0,1  . 60 

Fig.  40a  -  Exemple  de  niveaux  de  gris  avant  elimination  du  bruit 

Beispiel  von  Grauwerten  vor  der  Beseitigung  des  Rauschens . 61 

Fig.  40b  -  Niveaux  de  gris  apres  elimination  du  bruit  avec  un  nombre  seuil  egal  a  5 

Grauwerte  nach  der  Beseitigung  des  Rauschens  mit  einem  Schwellenwert  gleich  5 . 61 
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Fig.  41  -  Comparaison  des  erreurs  de  deplacement  avant  (courbes  reliant  les  croix) 

et  apres  elimination  du  bruit  (courbes  reliant  les  ronds)  ;  a)  courbes  obtenues 

en  autocorrelation  ;  b)  courbes  obtenues  en  correlation  croisee 

Vergleich  der  Verschiebungsfehler  vor  der  Beseitigung  des  Rauschens 

(Sternchen-Kurven)  und  nach  der  Beseitigung  des  Rauschens  (Kreis-Kun/en); 

a)  Kurven  bei  der  Autokorrelafion;  b)  Kun/en  bei  der  Kreuzkorrelation . 61 

Fig.  42  -  Pourcentages  des  simulations  donnant  des  erreurs  superieures  a  0,2  pixel 

pour  un  ecart  type  du  bruit  de  fond  egal  a  1  0  niveaux  de  gris  avant  elimination 

(courbes  du  haut)  et  apres  elimination  (courbes  du  bas)  du  bruit  ; 

a)  courbes  obtenues  en  autocorrelation  ;  b)  courbes  obtenues  en  correlation  croisee 

Prozentsafz  der  Simulationen,  die  bei  einer  Sfandardabweichung  des 

Rauschhintergrunds  gleich  7  0  Grauwerten  zu  Fehlern  von  uber  0,2  Pixel  fuhren, 

vor  der  Beseitigung  des  Rauschens  (Kurven  oben)  und  nach  der  Beseitigung 

des  Rauschens  (Kurven  unten);  a)  Kun/en  bei  der  Autokorrelafion;  b)  Kun/en  bei 

der  Kreuzkorrelation . . . 62 

Fig.  43  -  Pourcentages  des  simulations  donnant  des  erreurs  superieures  a  0,2  pixel 

pour  un  ecart  type  du  bruit  de  fond  egal  a  50  niveaux  de  gris  avant  elimination 
(courbes  du  haut)  et  apres  elimination  (courbes  du  bas)  du  bruit ;  a)  courbes  obtenues 
en  autocorrelation  ;  b)  courbes  obtenues  en  correlation  croisee 

Prozentsafz  der  Simulationen,  die  bei  einer  Sfandardabweichung  des  Rauschhintergrunds 
gleich  50  Grauwerten  zu  Fehlern  von  uber  0,2  Pixel  fuhren ,  vor  der  Beseitigung 
des  Rauschens  (Kurven  oben)  und  nach  der  Beseitigung  des  Rauschens  (Kun/en  unten); 
a)  Kurven  bei  der  Autokorrelafion;  b)  Kun/en  bei  der  Kreuzkorrelation . 62 

Fig.  44  -  Aire  d'interrogation  de  32  x  32  pixels  montrant  les  deplacements  imposes 

aux  particules  a  cause  du  gradient  d(Ax)/dy.  Le  vecteur  central,  a  I'ordonnee 
yo  =16  pixels,  re  presente  AXq 

Abfragebereich  von  32  x  32  Pixel,  der  die  Verschiebungen  der  Parfikel  aufgrund 

des  Gradienten  d(Ax)/dy  zeigt.  Der  Zentralvektor  bei  der  Ordinate  yo  —  16  Pixel 

sfellf  A Xq  dar . 64 

Fig.  45  -  Exemple  de  repartition  des  particules  correspondent  a  la  configuration  avec 

gradient  de  deplacement  d(AX)/dy  =  0,05  :  la  croix  correspond  a  la  position  des 
particules  a  I'instant  t]  et  le  rond  a  I'instant  t2 

Beispiel  der  Parfikelverteilung  entsprechend  der  Konfiguration  mif  Verschiebungs- 
gradienf  d(AX)/dy  =  0,05;  das  Sfernchen  entspricht  der  Position  der  Parfikel  zum 
Zeitpunkf  tj  und  der  Kreis  zum  Zeitpunkt  t2 . 64 

Fig.  46  a,b  -  Comparaison  des  histogrammes  suivant  I'axe  x  des  niveaux  de  gris  du  pic 

principal  de  la  fonction  d'autocorrelation  ;  a)  histogramme  obtenu  avec  presence 
d'un  gradient  de  deplacement ;  b)  histogramme  obtenu  sans  grandient  de  deplacement 
Vergleich  der  Saulendiagramme  gemaB  derx-Achse  der  Grauwerte  der  Hauptspitze 
der  Autokorrelationsfunktion:  a)  Saulendiagramm  bei  Vorhandensein  eines  Verschie- 
bungsgradienten;  b)  Saulendiagramm  ohne  Verschiebungsgradient . 65 

Fig.  46  c,d  -  Comparaison  des  histogrammes  suivant  I'axe  y  des  niveaux  de  gris  du  pic 

principal  de  la  fonction  d'autocorrelation  ;  c)  histogramme  obtenu  avec  presence  d'un 
gradient  de  deplacement ;  d)  histogramme  obtenu  sans  gradient  de  deplacement 
Vergleich  der  Saulendiagramme  gemaB  der  y-Achse  der  Grauwerte  der 
Fiauptspitze  der  Autokorrelationsfunktion;  c)  Saulendiagramm  bei  Vorhandensein  eines 
Verschiebungsgradienten;  d)  Saulendiagramm  ohne  Verschiebungsgradient . 65 

Fig.  47  -  Exemple  de  repartition  des  particules  correspondent  a  la  configuration  avec  gradient 
de  deplacement  d (AX) /d y  =  0,1 :  la  croix  correspond  a  la  position  des  particules  a 
I'instant  t]  et  le  rond  a  I'instant  t2 

Beispiel  der  Parfikelverteilung  entsprechend  der  Konfiguration  mif  Verschiebungs¬ 
gradient  d(AX)/dy  —  0,1 :  das  Sfernchen  entspricht  der  Position  der  Parfikel  zum 
Zeitpunkt  t j  und  der  Kreis  zum  Zeitpunkt  t2 . 66 
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Fig.  48  -  Comparaison  des  histogrammes  suivant  I'axe  x  des  niveaux  de  gris  du  pic 

principal  de  la  fonction  d'autocorrelation  ;  a)  histogramme  obtenu  avec  presence 

d'un  gradient  de  deplacement ;  b)  histogramme  obtenu  sans  gradient  de  deplacement 

Vergleich  der  Saulendiagramme  gemaB  derx-Achse  der  Grauwerte  der 

Haupfspitze  der  Autokorrelationsfunktion;  a)  Saulendiagramm  bei  Vorhandensein 

eines  Verschiebungsgradienten;  b)  Saulendiagramm  ohne  Verschiebungsgradienf . 67 

Legendes  tableaux  /  Tabellenlegenden 


Tab.  1  -  Caracteristiques  des  etudes  ISL  de  PIV 

Kennwerfe  der  Untersuchungen  des  ISL  miffels  PIV . .  1  7 

Tab.  2  -  Resultats  des  methodes  d'interpolation  sub-pixel 

Ergebnisse  der  Subpixel-lnterpolationsmefhoden . 37 

Tab.  3a  -  Erreur  de  de  deplacement  (moyenne  et  ecart  type)  en  pixel  pour  un  rayon  de  0,87  pixel 
Verschiebungsfehler  (mittlerer  und  Sfandardabweichung)  in  Pixel  bei  einem 
Radius  von  0,87  Pixel . 38 

Tab.  3b  -  Erreur  de  de  deplacement  (moyenne  et  ecart  type)  en  pixel  pour  un  rayon  de  1 ,60  pixel 
Verschiebungsfehler  (mittlerer  und  Sfandardabweichung)  in  Pixel  bei  einem 
Radius  von  1,60  Pixel . 38 
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d(AX)/dy  egal  a  0,05 
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Notations 


Lett  res  romaines 


Bj  j  niveau  de  bruit  de  Ny 

d]  distance  entre  plan  laser  et  objectif  de  la  camera  CCD 

62  distance  entre  objectif  et  matrice  CCD 

f  distance  focale  de  I'objectif 

f#  ouverture  relative  de  I'objectif 

G  gain  de  la  camera  CCD 

G(x,y,a)  fonction  gaussienne 

h(x,y)  reponse  impulsionnelle  en  amplitude 

(i/j)  coordonnees  du  pixel  courant 

1^  intensite  maximale  de  Pimage  de  la  particule  k 

I(x,y)  intensite 

M  grandissement  de  la  camera  CCD 

Nj  j  niveau  de  gris  du  pixel  (i,j) 

Rj<  rayon  de  la  particule  k 

R0,Ri/R.i  niveaux  de  gris  consideres  pour  I'interpolation  sub-pixel 
t  temps 

t]  instant  de  la  premiere  impulsion  laser 

t2  instant  de  la  seconde  impulsion  laser 

U (x/y/t)  amplitude  complexe 

V  vecteur  vitesse  dans  Pecoulement 


(x,y,z) 

JX/Y/Z) 

X 


coordonnees  dans  le  plan  image  (matrice  CCD) 
coordonnees  dans  le  plan  objet  (plan  laser) 
vecteur  position  dans  le  plan  objet 


Lettres  grecques 

5Z  profondeur  de  champ 

A  faille  d'un  pixel  de  matrice  CCD 

AT  intervalle  de  temps  entre  deux  impulsions  lasers  consecutives 
(Ax, Ay)  deplacement  impose  aux  particules 
Ax^  deplacement  entier  des  particules  suivant  I'axe  x 
Axe  deplacement  fractionnaire  des  particules  suivant  Paxe  x 
AX  vecteur  deplacement  d'une  particule  entre  les  instants  t  et  t+AT 
(j)  diametre  utile  de  I'objectif  photographique 

cp(kj)  fonction  d'autocorrelation  discrete 

X  longueur  d'onde  de  Pimpulsion  laser 

\|/(k,l)  fonction  de  correlation  croisee  discrete 
rayon  de  Pimage  de  la  particule  k 
V  frequence  spatiale 


Indices  et  abreviations 

cdm  centre  de  masse 

FFT  transformee  de  Fourier  rapide 

FTM  fonction  de  transfer!  de  modulation 

i  image  (matrice  CCD) 

g  geometrique 
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gau  gaussien 

h  reponse  impulsionnelle 

k  particule  k 

min  minimum 

max  maximum 

o  objet  (plan  laser) 

par  parabolique 

0  valeur  maximale 

*  operation  de  convolution 

~  proportionnalite 
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1 .  Introduction 


Depuis  quelques  annees,  la  technique  de  "Velocimetrie  par  Images  de  Par- 
ticules"  (ou  "Particle  Image  Velocimetry"  -  PIV  :  voir  [1,2])  s'affirme  comme 
une  technique  standard  de  mesure  de  champs  complets  de  vitesses  dans  les 
ecoulements.  A  I'lSL,  cette  technique  a  ete  mise  en  oeuvre  sur  diverses  ins¬ 
tallations  experimental  (tab.  1)  : 


Installation 

Theme  d'etude 

References 

bibliographiques 

soufflerie  subsonique 

decrochage  dynamique  sur  profil  oscillant 

[3]  a  [13] 

tunnel  hydrodynamique 

decrochage  dynamique  sur  profil  oscillant 

[14,15] 

soufflerie  supersonique 

diedre 

[14,15] 

soufflerie  supersonique 

culot  de  projectile 

[14,15] 

tunnel  de  tir  balistique 

projectile  en  vol 

[14,15] 

Tab.  1  -  Caracteristiques  des  etudes  I5L  de  PIV 


Pourtoutes  ces  etudes  (sauf  celle  concernant  le  decrochage  dynamique  sur 
profil  oscillant  en  soufflerie  subsonique  ou  une  methode  de  PIV  photogra- 
phique  avait  ete  utilisee),  les  mesures  ont  ete  effectuees  a  I'aide  d'un  systeme 
de  PIV  commercial  a  support  video  (enregistrement  sur  camera  CCD).  Afin 
d'utiliser  ce  systeme  dans  de  bonnes  conditions,  il  est  essentiel  d'avoir  une 
idee  de  la  precision  des  methodes  de  traitement  utilisees.  D'autre  part,  une 
analyse  de  la  litterature  scientifique  parue  sur  le  theme  de  la  PIV  montre  que 
ce  point  particulier  (la  precision  des  methodes  de  traitement  utilisees  en  PIV 
video)  n'a  ete  que  peu  etudie.  Enfin,  si  I'avantage  principal  de  ces  systemes 
commerciaux  de  PIV  est  de  permettre  la  realisation  de  mesures  dans  des 
conditions  "conviviales",  son  inconvenient  majeur  est  que  le  logiciel  de  de- 
pouillement  des  prises  de  vues  est  une  "boite  noire  fermee". 
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Introduction 


Ce  constat  fait  qu'il  est  essentiel  de  disposer,  a  I'lSL,  d'un  logiciel  de  de- 
pouiliement  de  prises  de  vues  "ouvert",  done  evolutif,  qui  serait  complemen- 
taire  du  logiciel  "ferme"  du  systeme  de  PIV  commercial  et  qui  permettrait 
d'etudier  la  precision  des  methodes  de  traitement  utilisees  en  PIV  video.  De 
plus,  ce  logiciel  pourrait  etre  ameliore  et  adapte  specifiquement  en  fonction 
de  I'ecoulement  etudie. 

* 

Pour  evaluer  les  performances  d'un  tel  logiciel  en  termes  de  precision  sur  la 
mesure  de  la  vitesse,  il  est  necessaire  de  disposer  d'un  modele  mathemati- 
que  permettant  de  simuler  des  images  de  PIV.  Un  tel  modele  a  ete  decrit  en 
detail  dans  [16].  Dans  le  present  Rapport,  apres  avoir  rappele  le  principe 
de  la  methode  de  PIV,  nous  exposerons  brievement  les  grandes  lignes  du 
modele  decrit  dans  [1  6],  puis  nous  I'appliquerons  a  diverses  configurations 
de  particules.  Les  simulations  numeriques  ainsi  effectuees  fournissent  des 
resultats  statistiques  importants  sur  la  precision  des  methodes  de  traitement 
utilisees  en  PIV  video. 
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2.  Moderation  d'un  systeme  de  PIV  video 


2.1  Rappel  du  principe  de  la  methode  de  PIV 

II  semble  que  la  methode  de  PIV  trouve  son  origine  dans  la  technique  du 
"speckle11  developpee  initialement,  vers  la  fin  des  annees  I960,  pour  des 
etudes  de  mecanique  des  solides.  Mais  le  vrai  depart  de  la  PIV  est  beau- 
coup  plus  recent  (fin  des  annees  1  980),  grace  principalement  au  fort  deve- 
loppement  de  la  mini-  et  micro-informatique.  De  nombreux  articles  de 
reference  sur  I'histoire  et  I'etat  de  Tart  de  la  PIV  pourront  etre  trouves  dans 

[1]  et  [13]. 

La  technique  de  PIV  reprend  la  majeure  partie  d'un  dispositif  classique  de 
visualisations  par  plan  laser  :  on  commence  done  par  creer  un  plan  de  lu- 
miere  laser  (fig.  1)  en  utilisant  des  optiques  appropriees  (lentilles  cylindri- 
ques  et/ou  spheriques).  On  ensemence  ensuite  I'ensemble  de  Pecoulement 
a  I'aide  de  petites  particules  (de  taille  typique  1  Jim)  et  on  enregistre  Pimage 
de  ces  particules  sur  un  support  photographique  ou  video.  Si  cette  opera¬ 
tion  est  repetee  un  court  instant  AT  (de  I'ordre  de  quelques  dizaines  de  jlis) 
plus  tard,  on  obtiendra  un  ensemble  de  doublets  de  deux  points,  chaque 
doublet  montrant  la  meme  particule  a  deux  instants  differents  (t  et  t  +  AT) 
et  deux  endroits  differents  (1(ft)  et  X(t+AT))  puisqu'elle  s'est  deplacee  entre 
les  prises  de  vues  en  suivant  Pecoulement  (fig.  2).  Si  I'on  parvient  ensuite  a 
regrouper  les  doublets  de  chaque  particule  et  a  trouver  leurs  positions  res- 
pectives,  on  determines,  a  partir  de  Pintervalle  de  temps  AT  separant  les 
deux  prises  de  vues,  la  norme  et  la  direction  du  vecteur  vitesse  local  par  : 

V(x,t)  MAT 

ou  M  represente  le  facteur  de  grandissement  entre  le  plan  image  (pellicule 
photographique  ou  matrice  CCD)  et  le  plan  objet  de  I'ecoulement. 
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Moc/e//safion  d’un  sysfeme  de  PIV  video 


Cette  operation  peut  etre  effectuee  automatiquement  par  une  technique  de 
correlation  locale  qui  va,  en  fait,  mesurer  le  deplacement  moyen  d'un  cer¬ 
tain  nombre  de  particules  se  trouvant  dans  une  "aire  d'interrogation"  (fig.  2) 
et  done  en  deduire  un  vecteur  vitesse  local  moyen  (voir  chap.  2). 


Fig .  2  -  Doublets  de  particules  dons  une  aire 
d'interrogation ;  (a)  en  blanc ;  les 
particules  correspondent  a  la  premiere 
impulsion  laser  et  en  noirf  a  la 
deuxieme ;  (b)  aspect  reel 


o  1st  image 
•  2nd  image 


II  reste  ensuite  a  determiner  le  signe  du  vecteur  vitesse  :  en  effet,  si  les  deux 
images  sont  enregistrees  sur  le  meme  support,  le  plus  souvent  photographi- 
que  (depouillement  par  autocorrelation),  il  n'est  pas  possible  de  savoir  le- 
quel  des  deux  points  correspond  au  premier  instant  de  prise  de  vues  et 
lequel  correspond  au  second  (fig.  2b  ci-dessus).  Pour  lever  cette  ambiguTte, 
comme  dans  le  cas  de  I'anemometrie  laser  Doppler  (LDA),  on  determine  la 
vitesse  par  rapport  a  une  reference  connue  non  nulle  :  ce  procede  est  connu 
sous  le  nom  de  decalage  d'image  ("image  shifting")  :  §  4.3.4  de  [1  3]. 

II  est  a  noter  ici  que  la  technique  de  PIV  peut  etre  employee  sans  avoir  re- 
cours  a  ('utilisation  d'un  dispositif  de  decalage  d'images.  En  effet,  I'ambiguY- 
te  sur  le  signe  de  la  vitesse  est  levee  des  lors  que  Ton  utilise,  par  exemple, 
deux  fenetres  de  correlation  differentes  provenant  de  deux  images  succes- 
sives  de  I'ecoulement  sur  une  camera  CCD  :  voir  [1  ].  Le  depouillement  des 
prises  de  vues  se  fait  alors  par  une  technique  d'intercorrelation  (ou  correla¬ 
tion  croisee)  qui  est  de  plus  en  plus  souvent  employee.  Dans  la  suite  de  ce 
Rapport,  nous  considererons  done  le  cas  de  prises  de  vues  par  camera 
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hAodellsoflon  dun  systeme  de  PIV  video 


CCD  (PIV  video)  et  les  deux  types  de  depouillement  possibles  :  autocorre¬ 
lation  et  correlation  croisee. 

Si  le  principe  general  de  la  methode  decrite  ci-dessus  semble  simple/  el le 
necessite  neanmoins  un  dispositif  experimental  tres  complet :  dispositif  d'en- 
semencement,  source  laser,  systeme  d'enregistrement  et  de  decalage  d'ima- 
ges.  L'ensemble  de  ce  dispositif  ainsi  que  certains  problemes  annexes 
(synchronisation,  choix  des  parametres  de  prises  de  vues,  validation  et  pre¬ 
cision  des  resultats,  etc.)  sonttraites  en  detail  dans  [1]  et  [1  3]. 

2.2  Modelisation  du  systeme  optique 

Comme  precise  dans  Pintroduction,  nous  ne  rappelons  ici  brievement  que 
les  grandes  lignes  du  modele  decrit  en  detail  dans  [1  6]. 

2.2.1  Notations 

La  figure  3  ci-dessous  schematise  de  maniere  simple  les  elements  princi- 
paux  dun  systeme  de  PIV  video  :  le  plan  objet  (plan  laser)  est  defini  par  un 
repere  direct  orthonorme  (0/X,Y,Z),  alors  que  le  plan  image  est  determine 
par  un  autre  repere  direct  orthonorme  (o,x,y,z)  tel  que  : 

•  o  et  O  se  trouvent  sur  la  normale  (axes  Z  et  z)  au  plan  de  Pobjectif  et 
sont  images  Pun  de  I'autre  par  cet  objectif, 

•  les  X  et  x  sont  paralleles  et  de  sens  oppose,  de  meme  que  les  axes  Y 
et  y, 

•  les  axes  Z  et  z  sont  confondus  et  de  meme  sens. 


plan  objet 
Objektebene 
plan  laser 
Lichtschnitt 


Fig.  3  -  Notation  pour  les  reperes  utilises 


Le  choix  de  ces  reperes  permet  d'eviter  les  signes  "moms"  qui  apparaissent 
d'habitude  dans  les  relations  analytiques  d'optique  geometrique  du  fait  des 
renversements  entre  objets  et  images  geometriques.  Avec  nos  notations, 
nous  aurons,  done  entre  objet  et  image  : 


x  =  MX  et 

ou  M  est  le  grandissement  image-objet  : 
En  utilisant  la  relation  de  conjugaison  : 


y  =  MY 


(2) 

(3) 

(4) 
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Modelisafion  d'un  systeme  de  PIV  video 


ou  f  est  la  distance  focale  de  I'objectif,  on  obtient  aussi  : 

(5)  d2  =  f(M+l)  e. 

Nous  supposerons,  pour  simplifies  I'intensite  I0(X,Y,Z)  du  plan  laser 
(IQ=  objet)  comme  etant  pratiquement  uniforme  et  constante  suivant  les 
axes  X  et  Y  (en  fait,  la  forme  de  la  loi  I0(X,Y)  ne  va  pas  intervenir  explicite- 
ment  dans  la  suite).  Suivant  I'axe  Z  (epaisseur  du  plan  laser),  nous  pren- 
drons  une  distribution  gaussienne  : 

(6)  I0(Z)  =  I0  exp  [-8(Z-Zq)2/AZq] 

ou  Z0  est  la  cote  du  centre  du  plan  laser  et  AZ0  son  epaisseur  correspon¬ 
dent  a  I0/eA  En  cela,  nous  nous  conformons  aux  modeles  adoptes  par 
d'autres  auteurs  [1  6]. 

Nous  supposerons,  d'autre  part,  que  I'epaisseur  AZ0  du  plan  laser  est  supe- 
rieure  a  la  profondeur  de  champ  du  systeme  optique  utilise,  donnee  par  : 

(7)  5Z  =  4(l+^)2f#2X  • 

Dans  cette  relation,  f#  designe  I'ouverture  relative  de  I'objectif  (f#  =  f/(|))  et 
X  la  longueur  d'onde  moyenne  du  faisceau  laser  creant  le  plan  laser.  Cette 
relation  est  obtenue  en  considerant  qu'un  objet  ponctuel  donne  une  image 
constitute  de  sa  tache  de  diffraction  theorique  (voir  la  demonstration  dans 
I'annexe  B  de  [13]).  A  Taide  de  cette  hypothese  AZ0  >  5Z,  nous  pourrons 
considerer  que  toutes  les  particules  a  I'interieur  du  plan  laser  donneront  des 
images  nettes  sur  la  matrice  CCD  et,  ainsi,  nous  n'aurons  pas  a  prendre  en 
compte  les  aberrations  de  mise  au  point. 

2.2.2  Relation  principale 

II  ressort  des  traites  d'optique  theorique  que,  de  fagon  generale,  I'image  fi¬ 
nale  d'un  objet  sur  un  capteur  optique  (pellicule  photographique, 
matrice  CCD)  est  donnee  par  la  convolution  entre  la  reponse  impul- 
sionnelle  de  la  pupille  de  sortie1  du  systeme  optique  et  I'image  geo- 
metrique  de  I'objet  (voir,  par  exemple  [1  7],  equation  6-5,  page  97). 

Cette  relation  generale  s'ecrit,  dans  le  cas  d'un  eclairement  coherent  en  ter- 
mes  de  champ  electrique  ou  magnetigne  U  reel  : 

(8)  Uj(x,y,t)  =  h(x,y)  *  Ug(x,y,t) 

ou  h(x,y),  t,  les  indices  i  et  g  designent  respectivement  la  reponse  impuision- 
nelle,  le  temps,  les  termes  d'image  et  de  geometrie.  L'intensite  correspon- 
dante  au  champ  ci-dessus  pourra  ainsi  se  calculer  suivant  : 


1 .  Pour  simplifier,  on  assimilera  la  pupille  de  sortie  a  la  forme  geometrique  de  I'objectif. 
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A/locfe/isaf/on  c/'un  s/sfeme  de  PIV  video 


Ij(x,y)  ~  |  h(x,y)  *  Ug(x,y,t)  1 2  (9) 

lei,  nous  n'explicitons  pas  les  differents  termes  multiplicatifs  constants  qui  in- 
terviennent  :  ainsi,  le  signe  —  designe  une  proportionnalite. 

En  PIV,  on  peut  raisonnablement  supposer  que  la  densite  de  particules  de 
Pecoulement  est  telle  que  Ion  peut  negliger  les  superpositions  des  taches 
centrales  de  diffraction  des  images  des  particules  dans  le  plan  image.  On 
peut  alors  montrer  que  la  relation  precedente  peut  s'ecrire  : 

ff+oo  2 

Ii(x,y)~  JJ  ^  I  h(x-xo,y-yo)  I  IgUo^O  dxgdyo  (10) 

=  I  h(x,y)  |  *  Ig(x,y)  • 

2.2.3  Modelisation  de  la  reponse  impulsionnelle  du  systeme  optique 

Dans  le  cas  d'un  objet  illumine  en  lumiere  monochromatique  coherente 
(cas  dun  laser),  la  reponse  impulsionnelle  d'une  lentil  le  mince  sans  aberra¬ 
tions  (done  uniquement  limitee  en  diffraction),  de  distance  focale  f  et  de  dia- 
metre  <[),  est  donnee  par  la  figure  de  diffraction  de  Frauenhofer  de 
I'ouverture  de  cette  lentille,  e'est-a-dire  la  figure  d'Airy  (voir  [1  7],  pp.  87-88) 
qui  s'exprime  a  I'aide  d'une  fonction  de  Bessel  de  premiere  espece  (voir 
[16]).  Une  approximation  particulierement  utile  pour  simplifier  les  calculs 
consiste  a  representer  la  figure  d'Airy  par  une  fonction  gaussienne  du  type  : 

G(x,y,o)  =  Gm  exp  [  ^  ^  ^ 

En  pratique,  dans  le  cas  d'un  objectif  reel  de  camera,  les  imperfections 
(aberrations)  doivent  etre  prises  en  compte.  Une  methode  simple  consiste  a 
quantifier  ces  aberrations  de  fagon  globale  par  I'utilisation  de  la  fonction  de 
transfert  de  I'optique  consideree  (voir,  par  exemple,  [1  7],  p.  105).  En  sup- 
posant  que  dans  le  cas  de  prises  de  vues  de  PIV  les  caracteristiques  optiques 
de  I'objectif  sont  bien  representees  par  la  FTM  (Fonction  de  Transfert  de 
Modulation)  mesuree  en  eclairement  incoherent  et  en  approximant  cette 
meme  FTM  par  une  fonction  gaussienne,  on  peut  finalement  montrer  que 
la  reponse  impulsionnelle  en  intensity  d'un  objectif  reel  sera  decrite  par  la 
fonction  gaussienne  G(x,  y,  a^)  ou  est  determine  par  la  FTM  de  I'objectif 
suivant  la  relation  : 


_ 1 _ 

°h  - 

ZTLOitm 


(12) 


Dans  cette  relation,  est  une  grandeur  caracteristique  de  la  FTM  de 
I'objectif  considere  (voir  [16]).  Si  Ton  suppose  enfin  que  la  reponse  impul¬ 
sionnelle  en  intensity  est  independante  de  la  position  de  la  particule  objet 
k  (invariance  spatiale  des  performances  optiques),  on  peut  alors  ecrire  la 
reponse  impulsionnelle  pour  une  particule  objet  situee  en  (X^Y^ZjJ  sous  la 
forme  : 
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(13)  |  h(x,y)  I  2  =  |  ho  1 2  GCx-MX^y-MY k.ah)  avec  Oh  =  - • 

ZJlGpTM 

2.2.4  Moderation  de  I'image  geometrique  de  I'objet 

Dans  le  cas  d'un  eclairement  coherent,  I'amplitude  complexe  de  I'image 
geometrique  d'un  objet  est  defini  dans  le  plan  image  (x,y)  par  (voir  [17], 
equ.  5-43,  p.  89)  : 

0  4)  Ug(x,y)  =  mUo(m  ’  M-* 

ou  U0(X,Y)  represente  I'amplitude  de  I'onde  complexe  immediatement  der- 
riere  I'objet  (voir  [1  7],  p.  84). 

En  supposant,  comme  au  paragraphe  2.2.2,  que  les  taches  centrales  de 
diffraction  de  ces  particules  ne  se  superposent  pas  dans  le  plan  image,  I'in- 
tensite  dans  le  plan  image  s'ecrit,  pour  un  objet  constitue  de  n  particules  : 


(15) 


-Lvl  TT  .  X,  |2 


Ig(x,y) = X '  uo,k(; 


k=l 


M’M 


— Yi  r—  -V 
m2  Z  ’  W 

1V1  k=l 


Pour  des  raisons  de  simplicity  dans  I'operation  de  convolution  (la  convolu¬ 
tion  de  deux  fonctions  gaussiennes  reste  une  fonction  gaussienne)  et  puis- 
que  la  reponse  impulsionnelle  de  I'optique  est  elle-meme  de  forme 
gaussienne,  nous  representerons  egalement  chaque  terme  IQ  ^  par  des 
fonctions  gaussiennes,  et  I'intensite  de  I'image  geometrique  s'ecrira  finale- 
ment  : 

n 

1  V1  MRk 

(1 6)  Ig(x,y)  =  ~ ”2^  %  kt  G(x-MXk  ,  y-MYk  ,  — )  * 

k=l 


Cette  formule  permet  ainsi  de  calculer  I'image  geometrique  dans  le  plan 
image  (o,x,y)  d'un  ensemble  de  n  particules  de  rayon  situes  aux  points 
(Xk,Yk,Zk). 

2.2.5  Calcul  de  I'  intensity  lumineuse  arrivant  sur  la  matrice  CCD 

Dans  le  cas  d'un  objectif  ideal  uniquement  limite  en  diffraction,  nous  pou- 
vons  ecrire,  d'apres  les  paragraphes  precedents,  ('intensity  lumineuse  arri¬ 
vant  sur  la  camera  CCD  sous  la  forme  : 

(1 7)  Ii(x,y)  =  I  h(x,y)  | 2  *  Ig(x,y) 

ou  le  terme  |h(x,y)P  est  represente  par  une  fonction  gaussienne.  Or,  la 
convolution  de  deux  fonctions  gaussiennes  reste  une  fonction  gaussienne 
et,  par  consequent,  nous  obtenons  finalement  : 


< 
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n 


V  Ck 

Ii(x,y)  =  2^1  k  G(x-MXk  ,  y-MYk  ,  — ) 

(18) 

k=l 

avec 

<?k  =^/(MRk)2  +  (l,22Xf#(M+l))2 

(19) 

et 

tc*2  Rk  t 

Ik=  4  I0,k  kt(-)2- 

(20) 

Dans  cette  expression,  le  terme  1^  represente  I'intensite  maximale  de  I'ima- 
ge  de  la  particule  k  et  le  terme  le  rayon  en  intensite  de  I'image  de  cette 
particule  (valeur  a  \/e  ). 


Dans  le  cas  d'un  objectif  reel,  nous  avons  de  la  meme  fagon  que  ci-dessus  : 

n 

Ii(x,y)  =  ^  Ik  G(x-MXk  ,  y-MYk  ,  y)  (21) 

k=l 


mais  avec  cette  fois: 


et 


ak  =  A  / (MRk)2  +  C 


1 


L  =■ 


JcahTM 

Rk  - 


r 


2kg, 


2  l0,k  kt  (Ck) 


FTM 


(22) 

(23) 


Comme  au  paragraphe  precedent,  le  terme  1^  represente  I'intensite  maxi¬ 
male  de  I'image  de  la  particule  k  et  le  terme  G^,  le  rayon  en  intensite  de 
I'image  de  la  particule  k  (valeur  a  \/e). 

2.3  Moderation  de  la  camera  CCD 

Nous  avons  vu  au  paragraphe  precedent  que,  dans  tous  les  cas,  I'intensite 
lumineuse  diffusee  par  les  particules  et  arrivant  sur  la  matrice  CCD  est  don- 
nee  par  la  relation  suivante  : 

n 

Ii(x,y)  =  ^  Ik  G(x-MXk  ,  y-MYk  ,  y)  (24) 

k=l 


ou  les  expressions  de  et  varient  suivant  I'approximation  consideree 
(objectif  ideal  ou  reel).  II  s'agit  maintenant  de  modeliser  I'action  de  la  ma¬ 
trice  CCD  pour  en  deduire  le  signal  delivre  a  sa  sortie.  Nous  nous  limiterons 
aux  caracteristiques  principales  d'une  matrice  CCD  schematisees  par  les 
differents  points  suivants  : 
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•  integration  spatiale  sur  chaque  pixel, 

•  gain  et  bruit, 

•  quantification. 

Precisons  ici  la  disposition  de  la  matrice  CCD  :  nous  supposons  celle-ci 
constitute  de  n  X  n  pixels  selon  le  schema  de  la  figure  4  (le  point  de  coor- 
donnees  (i,j)  correspond  au  coin  inferieur  gauche  du  pixel  (i,j)). 

n 

j 


Fig.  4  -  Disposition  et  notation  pour  les  pixels  de  la 
matrice  CCD 


1 

1  i  n 


2.3.1  Integration  spatiale 

Une  des  caracteristiques  principales  d'une  camera  CCD  est  de  donner, 
pour  chaque  pixel  (i,j),  un  signal  de  sortie  proportionnel  a  Tintensite  moyen- 
ne  re^ue  par  ce  pixel.  Ceci  peut  s'ecrire  de  faqon  simplifiee  sous  la  forme 
suivante  : 

(25)  sij~77tt  Ii(x,y)dxdy 

(i.j) 

ou  A  designe  la  taille  d'un  pixel  (suppose  etre  carre).  Avec  I'expression  ci- 
dessus  de  Ij(x,y),  il  vient  finalement  : 

(26)  Sy  -  f  (p2  I  Ik  [erf(T2)  -  erf(T, )]  [erf(T4)  -  erf(T3 )]  • 

k=1 


Les  integrales  qui  interviennent  ci-dessus  peuvent  s'exprimer  aisement  par 
I'internnediaire  de  la  "fonction  d'erreur"  erf(x)  definie  par : 


[0+l)A  -  MXk]  ; 
[(j+I)A-MYk]. 


2.3.2  Gain  et  bruit 

Pour  simplifier,  considerons  que  le  seul  reglage  de  la  camera  CCD  est  celui 
d'un  gain  G.  De  plus,  le  bruit  introduit  par  les  circuits  de  la  camera  sera 
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note  Bj  j  our  chaque  pixel  (i,j)  de  telle  sorte  que  nous  ecrirons  le  signal  de 
sortie  pour  ce  pixel  (i,j)  sous  la  forme  : 

sg=GsiJ+Bij'  '29> 

2.3.3  Quantification 

Dans  le  cas  d'une  camera  CCD  classique,  la  quantification  s'effectue  ge- 
neralement  sur  8  bits,  c'est-a-dire  sur  256  niveaux  de  gris  compris  entre  0 
et  255.  En  assimilant  le  signal  d'entree  du  quantificateur  a  S7^,  le  niveau 
de  gris  Nj  j  du  pixel  (i,j)  s'ecrira  finalement  : 


(30) 


Nj ;  =  0  si  S. .  <0 

1J  i.J 

i  i 

Njj  =  E(  S.  .  +  1/2)  si  0  <  S.  <  255 

t 

Nil  =255  siS. .  >255 
1J  iJ 


Le  cas  ideal  est  celui  ou  le  bruit  Bj  j  est  nul  ou  negligeable  et  ou  le  gain  est 
parfaitement  regie  de  telle  sorte  que  la  dynamique  complete  du  quantifica¬ 
teur  8  bits  soit  utilisee,  soit  N^min  =  0  et  (Nj  j)max  =  255.  En  fait,  les  prises 
de  vues  reelles  de  PIV  sont  souvent  loin  de  satisfaire  a  ces  conditions  (voir 
[16]). 


2.4  Validation  du  modele  analytique 

Le  modele  analytique  developpe  ici  necessite  d’etre  valide  avantson  utilisa¬ 
tion  pour  des  simulations  numeriques.  On  peut  alors  faire  les  commentaires 
suivants  : 

•  les  autres  modeles  de  PIV  proposes  jusque-la  font  tous  I'hypothese 
d'images  de  particules  dont  I'intensite  est  de  repartition  gaussienne 
(voir  [1,2]),  comme  dans  notre  modele  ; 

•  pour  la  quasi-totalite  des  autres  modeles  de  PIV  proposes  jusque-la, 
on  suppose  le  systeme  optique  de  prises  de  vues  comme  etant  parfait, 
sans  aberrations,  ce  qui  mene  a  la  relation  (19)  pour  le  rayon  des 
images  de  particules.  Dans  [1  6],  on  montre  que  cette  relation  sous- 
estime  le  rayon  des  images  de  particules,  alors  que  la  relation  (22) 
en  donne  une  estimation  correcte. 
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Comme  pour  le  chapitre  precedent,  seules  les  relations  principales  sont 
rappelees  ici  :  les  details  sont  explicites  dans  [1  6]. 

3.1  Autocorrelation  et  correlation  croisee 

Dans  ce  paragraphe,  nous  n'abordons  que  les  formes  les  plus  simples  des 
methodes  de  correlation  utilisees  en  PIV.  D'autres  formes  de  correlation 
existent  et  sont  classiquement  utilisees  en  PIV  (par  exemple,  le  coefficient  de 
correlation  :  voir  [1  6]). 

3.1.1  Autocorrelation 

Classiquement,  le  traitement  des  prises  de  vues  de  PIV  est  effectue  par  une 
technique  d'autocorrelation  lorsque  les  deux  prises  de  vues  aux  instants  t  et 
At+T  sont  enregistrees  sur  la  meme  image  video.  On  calcuie  alors  une 
fonction  discrete  d'autocorrelation  qui  peut  s'ecrire  sous  sa  forme  la  plus 
simple  : 

cp(k,  1)  =  £  £  Ni(j  Ni  +  k  j  +  1  (3i) 

i  j 

ou  les  Ni,j  sont  calcules  a  partir  du  modele  du  chapitre  2.  La  figure  d'auto¬ 
correlation  ainsi  obtenue  sera  composee  d'un  pic  central  tres  eleve  corres¬ 
pondent  a  un  decalage  k  =  I  =  0  (pic  d'autocorrelation  d'ordre  0)  et  deux 
pics  (d'ordre  1)  symetriques  par  rapport  a  ce  pic  central  et  dont  la  position 
par  rapport  au  pic  central  correspond  a  la  moyenne  des  deplacements  de 
toutes  les  particules  de  I'aire  d'interrogation.  La  vitesse  locale  moyenne  de 
I'ecoulement  dans  cette  aire  s'obtient  ensuite  aisement  en  determinant  nu- 
meriquement  la  position  de  ces  pics  et  en  tenant  compte  de  I'intervalle  de 
temps  AT  entre  les  deux  impulsions  lasers  et  du  grandissement  M  entre  le 
plan  objet  et  la  matrice  CCD,  selon  ('equation  (1). 

Les  termes  (p(k,I)  peuvent  etre  calcules  directement  selon  I'equation  (31), 
mais  il  est  bien  plus  rapide,  d'un  point  de  vue  numerique,  d'utiliser  des 
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transformees  de  Fourier  rapides  (FFT)  directes  et  inverses  (voir  le  theoreme 
de  Wold  et  Wiener-Khintchine  [16]). 

3.1 .2  Correlation  croisee 

Lorsque  les  deux  prises  de  vues  aux  instants  t  et  t  +  AT  sont  enregistrees  sur 
deux  images  video  differentes,  le  depouillement  se  fait  classiquement  par 
une  technique  de  correlation  croisee  et  on  calcule  alors  une  fonction  dis¬ 
crete  de  correlation  croisee  qui  peut  s'ecrire  sous  sa  forme  la  plus  simple  : 

(32)  ¥(k.l)=s2INifjNi+k,j  +  I 

i  j 

ou  Nj  j  designe  le  niveau  de  gris  de  I'image  enregistree  a  I'instant  t  et  N'jj, 
le  niveau  de  gris  de  I'image  enregistree  a  I'instant  t  -f  AT.  La  figure  de  cor¬ 
relation  croisee  ainsi  obtenue  sera,  cette  fois,  composee  d'un  seul  pic  eleve 
dont  la  position  par  rapport  au  centre  de  la  figure  de  correlation  correspond 
a  la  moyenne  des  deplacements  de  toutes  les  particules  de  I'aire  d'interro- 
gation.  La  vitesse  locale  moyenne  de  I'ecoulement  dans  cette  aire  s'obtient 
ensuite  de  la  meme  fagon  que  dans  le  cas  de  ('autocorrelation. 

La  encore,  les  termes  \|/(k,I)  peuvent  etre  calcules  directement  selon  ['equa¬ 
tion  (32),  mais  il  est  plus  rapide,  d'un  point  de  vue  numerique,  d'utiliser  des 
transformees  de  Fourier  rapides  (FFT)  directes  et  inverses. 

3.2  Interpolation  sub-pixel 

Les  methodes  d'autocorrelation  et  de  correlation  croisee  permettent  de  re- 
trouver  le  nombre  entier  de  pixels  (appele  "deplacement  entier")  le  plus  pro- 
che  des  vrais  deplacements  programmes  (la  precision  de  mesure  est  done 
de  ±  1/2  pixel).  La  precision  de  mesure  d'un  systeme  de  PIV  video  sera 
alors,  dans  une  certaine  mesure,  independante  de  ce  deplacement  entier  et 
ne  dependra  que  du  deplacement  fractionnaire  impose. 

Pour  obtenir  une  meilleure  precision,  il  est  necessaire  de  trouver  le  "depla¬ 
cement  fractionnaire"  manquant  :  des  methodes  dites  "d'interpolation  sub- 
pixel"  sont  necessaires.  Pour  la  composante  de  deplacement  Ax  suivant 
I'axe  x,  nous  ecrivons  : 

(33)  Ax  =  AxE  +  Axe 

ou  Axg  et  Axe  designent  respectivement  les  deplacements  "entier"  et  "frac- 
tionnaire".  Dans  la  suite,  nous  nous  attacherons  a  determiner  Axe. 

II  existe  actuellement,  dans  la  litterature,  trois  methodes  classiques  pour 
trouver  Axe  : 

•  une  methode  de  centre  de  masse  (barycentre), 

•  une  methode  d'ajustement  de  fonction  parabolique, 

•  une  methode  d'ajustement  de  fonction  gaussienne. 
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Pour  ces  trois  methodes,  la  premiere  etape  consiste  a  effectuer  un  change- 
ment  de  repere  par  translation  sur  les  axes  x  et  y  de  la  figure  de  correlation 
consideree  de  telle  sorte  que  le  pixel  de  niveau  de  gris  le  plus  eleve  (note 
Rq)  se  trouve  a  I'origine  x  =  y  =  0  de  ce  nouveau  repere.  Nous  ne  consi- 
dererons  ensuite  que  les  niveaux  de  gris  immediatement  a  gauche  et  a  droi- 
te  de  Rq,  niveaux  qui  seront  respectivement  notes  R_]  et  R]  (fig.  5).  Les  autres 
niveaux  de  gris  sont  consideres  comme  immerges  dans  le  "bruit  de  fond"  de 
la  figure  de  correlation  et  n'interviendront  done  plus  dans  la  suite.  Le  depla¬ 
cement  fractionnaire  Axe  sera  alors  donne  par  I'abscisse  du  "centre"  de  la 
distribution  des  niveaux  de  gris  Rq  ^  Rq  et  R] .  Les  trois  methodes  developpees 
dans  la  suite  se  distinguent  par  la  fagon  de  definir  le  "centre"  de  ces  trois 
niveaux  de  gris.  On  obtient  finalement  (voir  [16])  : 


•  Methode  de  centre  de  masse  (barycentre)  : 

Ri-R-i  . 
(Ax£)cdm  “  R_j  +  R0  -f  Rj 

•  Methode  d'ajustement  de  fonction  parabolique  : 

R-i-Ri 

(Axe)par  “  2(R  ,  +  Rj  -  2R0) 

•  Methode  d'ajustement  de  fonction  gaussienne  : 

In  R_i  -  In  R| 

(Axe)gau  ~  2(ln  R_,  +  In  R,  -  21n  R0) 


(34) 


(35) 


(36) 


Fig.  5  -  Methode  d'ajustement  par  des 
fonctions  paiabolique  et 
gaussienne 


axe  x 
x-Achse 
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4.1  Generalites 

Dans  ce  Rapport,  plusieurs  configurations  de  particules  seront  considerees  : 

•  dans  un  premier  temps,  nous  simulerons  des  prises  de  vues  sans  bruit 
de  fond  :  le  premier  des  deux  paragraphes  concernes  traitera  le  cas 
d'une  seule  particule  presente  dans  I'aire  d'interrogation  (§  4.2),  tan- 
dis  que  le  second  traitera  le  cas  de  plusieurs  particules  presentes  dans 
cette  aire  d'interrogation  (§  4.3)  ; 

•  dans  une  seconde  partie,  nous  introduirons  du  bruit  dans  les  simula¬ 
tions  des  differentes  prises  de  vues  (§  4.4)  ; 

•  enfin,  dans  une  troisieme  partie,  nous  aborderons  le  cas  d'a ires  d'in¬ 
terrogation  en  presence  de  gradients  de  deplacement  (§  4.5). 

Avant  de  commencer  les  simulations,  il  est  necessaire  de  preciser  quelques 
points. 

4.1.1  Variables  aleatoires 

La  methode  de  PIV  n'est  pas  une  methode  entierement  determinate  dans  le 
sens  ou  elle  comporte  un  caractere  aleatoire  important  du  principalement 
aux  particules  utilisees.  En  effet,  en  pratique,  les  grandeurs  liees  aux  parti¬ 
cules  (position,  taille,  densite,  deplacement)  sont  des  grandeurs  aleatoires 
et  ne  peuvent  etre  definies  que  de  fagon  statistique.  Precisons  ces  differents 
points  : 

•  position  des  particules  :  elle  depend  du  generateur  de  particules  et 
de  I'ecoulement.  En  toute  rigueur,  la  position  initiale  des  particules 
(i.e.  lors  de  la  premiere  exposition)  est  fonction,  pour  un  generateur 
de  particules  donne,  de  la  turbulence  de  I'ecoulement  considere. 
Pour  simplifier,  nous  representerons  la  position  initiale  des  particules 
par  une  distribution  uniforme  ; 
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•  rayon  des  particules  :  le  rayon  reel  des  particules  depend  du  genera¬ 
tes  de  particules  et  est  souvent  represente  par  une  distribution  log- 
normale  ou  de  Weibull  [13].  Dans  nos  simulations,  ce  n'est  pas  le 
rayon  reel  des  particules  qui  intervient,  mais  celui  des  images 
des  particules.  Pour  simplifier,  nous  representerons  le  rayon  par 
une  distribution  normale  (gaussienne)  ; 

•  densite  de  particules  :  elle  depend  egalement  du  generateur  de  par¬ 
ticules  et  de  I'ecoulement  et  sera  representee  par  une  distribution  nor¬ 
male  ; 

•  deplacement  des  particules  :  il  depend  de  I'ecoulement  et  sera  (pour 
une  aire  d'interrogation  donnee)  soit  constant  pour  un  ecoulement 
uniforme  (§  4.2  a  4.4),  soit  variable  dans  le  cas  de  presence  de  gra¬ 
dients  de  deplacement  (§  4.5). 

Une  autre  grandeur  aleatoire  est  constitute  par  I'ensemble  des  sources  de 
bruit  du  dispositif  de  prises  de  vues  (camera  CCD,  reflexions  parasites,  etc.) 
et  sera  modelisee  au  paragraphe  4.4  par  une  distribution  normale. 

4.1.2  Conditions  optimales  en  PIV 

Des  etudes  anterieures  ont  determine  les  conditions  optimales  de  prises  de 
vues  de  PIV  [1 ,2,1  3].  Ces  conditions  portent  essentiellement  sur  le  nombre 
de  particules  presentes  dans  chaque  aire  d'interrogation,  sur  leur  deplace¬ 
ment  et  leur  rayon  : 

•  le  nombre  de  particules  presentes  dans  la  fenetre  d'interrogation  (de 
taille  32  pixels  x  32  pixels,  dans  notre  cas)  doit  etre  superieur  a  7  en 
correlation  croisee,  et  de  I'ordre  de  1  5  en  autocorrelation  ; 

•  le  deplacement  des  particules  doit  etre  de  I'ordre  du  quart  ou  du  tiers 
de  la  taille  de  la  fenetre  d'interrogation  en  autocorrelation  (dans  no¬ 
tre  cas,  compris  entre  8  et  1  0  pixels). 

4.1 .3  Modele  numerique 

Une  synthese  des  relations  analytiques  obtenues  au  chapitre  2  montre  que 
la  repartition  des  niveaux  de  gris  de  chaque  pixel  depend  directement  des 
trois  variables  suivantes  (la  taille  des  pixels  de  la  CCD  etant  une  constante)  : 

•  les  positions  x^  =  MXj<  et  y^  =  MY^  des  centres  des  images  de  parti¬ 
cules  sur  la  CCD  :  elles  seront  representees  par  une  distribution  uni¬ 
forme  (voir  §  4.1 .1)  ; 

•  le  rayon  des  images  de  ces  particules  :  pour  des  prises  de  vues  ty- 
piques  de  PIV,  on  obtient  des  rayons-images  de  I'ordre  d'une  dizaine 
de  |Llm,  ce  qui  correspond  a  1  a  2  pixels  ; 

•  I'intensite  maximale  de  I'image  de  ces  particules,  determinee  par  le 
produit  G  G  etant  ici  le  gain  de  la  CCD.  Ce  produit  G  1^  fait  in¬ 
tervene  un  certain  nombre  de  grandeurs  inconnues  (Ig  et  kt).  Pour 
simplifier,  nous  considererons  le  terme  G  \  comme  constant. 

A  partir  de  ce  qui  precede,  nous  sommes  maintenant  en  mesure  de  cons- 
truire  un  modele  numerique  pour  simuler  un  systeme  de  PIV  video. 
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4.2  C as  cTune  seule  particule  presente  dans  Voire 
d'interrogatio  n 

Pour  commencer,  considerons  le  cas  le  plus  simple  possible  :  a  1‘instant  t-j 
de  la  premiere  prise  de  vue  de  PIV  ne  se  trouve  qu'une  seule  particule  dans 
I'aire  d'interrogation.  A  I'instant  +2  —  t]  +  AT,  cette  meme  particule  se  trouve 
toujours  dans  la  meme  aire  d'interrogation,  mais  dans  une  autre  position. 
Nous  supposons  qu’aucune  autre  particule  n'entre  dans  cette  aire  d'interro¬ 
gation  entre  les  instants  t]  et  t2-  D'autre  part,  nous  supposons  que  les  con¬ 
ditions  de  prises  de  vues  et  de  reglage  de  la  camera  CCD  sont  parfaites, 
c'est-a-dire  (voir  §  4.1)  : 

•  pas  de  bruit  :  B^j  =  0, 

•  CCD  juste  saturee  :  G  =  255. 

4.2.1  Exemple 

La  figure  6  montre  un  exemple  de  repartition  tridimensionnelle  des  niveaux 
de  gris  N{  j  dans  le  cas  considere  ici.  Cette  figure  a  ete  calculee  a  partir  des 
relations  (30)  du  paragraphe  2.3.3  et  represente  la  meme  particule  aux  ins¬ 
tants  t]  et  t2-  Les  niveaux  de  gris  de  la  figure  6  ont  ete  adimensionnes  par 
255. 


Les  positions  (x,y)  et  les  deplacements  (Ax, Ay)  de  la  particule  (choisis  au  ha- 
sard)  sont  : 

•  x(t-| )  =  10,21  pixels  ety(t])  =  12,86  pixels, 

•  x(t2)  —  22,92  pixels  et  y(t2)  =  1  7,00  pixels, 

•  Ax  =  x(t2)  -  x(t] )  =  12,71  pixels  et  Ay  =  yfo)  -  y(t-j )  =  4,14  pixels. 

Le  rayon  de  la  particule  consideree  a  ete  fixe  a  C]  =  0,87  pixel.  On  remar- 
quera  que  le  deplacement  de  la  particule  est  superieur  au  deplacement  op¬ 
timal  indique  au  paragraphe  4.1.2,  mais  ceci  n'a  ici  aucune  importance 
puisqu'une  seule  particule  intervient  ici. 

La  figure  7  represente  la  fonction  d'autocorrelation  correspondant  a  la  figu¬ 
re  6.  On  obsen/e  bien  un  pic  central  eleve  et  deux  pics  principaux  symetri- 
ques  par  rapport  a  ce  pic  central,  et  qui  donne  le  deplacement  des 
particules  [1  6]. 
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Fig .  7  -  Fonction  d't autocorrelation  correspondent  6  lo 
figure  6 


Dans  le  cas  ou  les  deux  prises  de  vues  a  t]  et  ^  sont  enregistrees  sur  deux 
images  differentes  (depouillement  par  correlation  croisee),  on  obtient  les  fi¬ 
gures  8a  et  8b. 


Fig.  8a  -  Repartition  tridimensionnelle  des  niveaux  de  gris  d  I'intant  t]  Fig.  8b  -  Repartition  tridimensionnelle  des  niveaux  de  gris  a  I'instant  t2 


La  fonction  de  correlation  croisee  correspondante  est  representee  sur  la  fi¬ 
gure  9.  ['application  des  methodes  d'interpolation  sub-pixel  a  cette  confi¬ 
guration  de  particules  permet  alors  de  determiner  precisement  le 
deplacement  des  particules  [1  6].  Le  tableau  2  resume  alors  les  resultats  ob- 
tenus.  Dans  ce  tableau,  les  colonnes  Ax  et  Ay  designent  respectivement  les 
deplacements,  suivant  les  axes  x  et  y,  obtenus  par  les  methodes  de  correla¬ 
tion  et  d'interpolation  sub-pixel  alors  que  les  colonnes  ex  et  ey  indiquent  res¬ 
pectivement  les  differences  : 

ex  =  Ax  -  1  2,71  (en  pixel), 
ey  =  Ay  -  4,1  4  (en  pixel), 

soit  les  differences  entre  les  deplacements  obtenus  par  les  methodes  de  cor¬ 
relation  et  les  deplacements  imposes.  Ces  differences  seront,  dans  la  suite, 
appelees  “erreurs  de  deplacement"  et  seront  la  principale  grandeur  etudiee 
dans  la  suite  de  ce  Rapport.  Notons  que,  dans  ce  cas,  d'une  seule  particule 
par  aire  d'interrogation,  sans  bruit,  les  resultats  obtenus  ne  dependent  pas 
de  la  methode  de  correlation,  mais  uniquement  de  la  methode  d'interpola¬ 
tion  sub-pixel  choisie. 


36 


ISL  -  R  119/98 


S/mu/afions  numer/ques 


Ax  [pixel] 

ex  [p>xel] 

Ay  [pixel] 

ey  [pixel] 

centre  de  masse 

12,74 

+0,03 

4,13 

+0,01 

parabole 

12,67 

-0,04 

4,14 

0,00 

gaussienne 

12,61 

-0,10 

4,19 

-0,05 

Tab.  2  -  Resultats  des  methodes  d'interpolation  sub-pixel 


4.2.2  Convergence  de  I'erreur  de  deplacement 

L'erreur  de  deplacement,  definie  a  la  fin  du  paragraphe  precedent,  sera 
maintenant  etudiee  de  maniere  statistique.  Pour  cela,  tragons  les  courbes 
de  convergence  de  I'erreur  de  deplacement  (en  moyenne  et  en  ecart  type) 
en  fonction  du  nombre  de  simulations  (300,  en  tout)  pour  une  particule  de 
rayon  egal  a  0,87  pixel  dont  la  position  a  Tinstant  t]  est  aleatoire  et  le  de- 
placement  identique  a  celui  de  I'exemple  du  paragraphe  4.2.1  (fig.  1  0a  et 
10b). 


Zahl  von  Iterationen  Zahl  von  Iterationen 


Fig .  10a  -  Convergence  de  I'erreur  de  deplacement  ( moyenne ,  en 
pixel)  en  fonction  du  nombre  de  simulations  pour  les  trois 
methodes  d'interpolation  sub-pixel  (rayon  de  0,87  pixel) 


Fig.  10b  -  Convergence  de  I'erreur  de  deplacement  (ecart  type ,  en 
pixel)  en  fonction  du  nombre  de  simulations  pour  les  trois 
methodes  d'interpolation  sub-pixel  (rayon  de  0,87  pixel) 
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Les  convergences  de  I'erreur  en  moyenne  et  en  ecart  type  sont  etablies  assez 
rapidement.  II  suffit  d'effectuer  une  centaine  d'iterations  pour  obtenir  des 
resultats  satisfaisants  :  au-dela  de  cette  valeur,  ni  I'erreur  moyenne,  ni 
I'ecart  type  ne  varient  de  plus  d'un  centieme  de  pixel. 

Tragons  les  memes  courbes  pour  un  rayon  de  1 ,6  pixel  (fig.  1  1  a  et  1 1  b).  La 
convergence  est  obtenue  beaucoup  plus  rapidement  que  pour  le  rayon  de 
0,87  pixel.  Par  contre,  la  valeur  moyenne  de  I'erreur  obtenue  par  la  metho- 
de  d'interpolation  par  le  centre  des  masses  augmente  beaucoup  lorsque  le 
rayon  passe  de  0,87  a  1,6  pixel. 
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0,16 

0,14 

0,12 

0,1 
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0,06 

0,04 
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o  o 


moyenne  sur  300  iterations 
Mittelwert  aus  300  Iterationen 
o :  cdm 
+ :  par 
x :  gau 
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Fig.  77a- Convergence  de  I'erreur  de  deplacement  ( moyenne ,  en 
pixel)  en  fonction  du  nombre  de  simulations  pour  les  trois 
metbodes  d'interpolation  sub-pixel  (rayon  de  1,60  pixel) 


Fig.  lib-  Convergence  de  I'erreur  de  deplacement  (ecart  type,  en 
pixel)  en  fonction  du  nombre  de  simulations  pour  les  trois 
metbodes  d'interpolation  sub-pixel  (rayon  de  1,60  pixel) 


Les  resultats  ainsi  obtenus  sont  consignes  dans  les  tableaux  3a  et  3b  sui- 
vants  : 


Methode 

Valeur  moyenne 

Ecart  type 

centre  des  masses 

0,056 

0,020 

parabolique 

0,066 

0,072 

gaussienne 

-0,013 

0,058 

Tab.  3a  -  Frreur  de  deplacement  (moyenne  et  ecart  type)  en  pixel  pour  un  rayon  de  0,87  pixel 


Methode 

Valeur  moyenne 

Ecart  type 

centre  des  masses 

0,172 

0,001 

parabolique 

0,033 

1 0,003 

gaussienne 

0,000 

0,003 

Tab.  3b  -  Erreur  de  deplacement  (moyenne  et  ecart  type)  en  pixel  pour  un  rayon  de  1,60  pixel 


4.2.3  Repartition  de  I'erreur  de  deplacement 

Pour  tracer  les  courbes  de  convergence  precedentes,  nous  avons  effectue 
300  simulations.  Reprenons  ces  memes  resultats  pour  tracer  les  histogram- 
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mes  representant  la  repartition  de  I'erreurde  deplacement  commise  suivant 
Taxe  x  (en  nombres  de  pixels)  lorsque  la  particule  presente  a  un  rayon  de 
0,87  pixel  (fig.  1  2a  a  1 2c).  Chaque  figure  correspond  a  une  methode  d'in- 
terpolation  parmi  les  trois  proposees  (centre  des  masses,  parabole,  gaus- 
sienne). 


a)  methode  d'interpolation :  centre  des  mosses  b)  methode  d'interpolation :  porobole 


c)  methode  d'interpolation :  goussienne 


Fig.  12  -  Histogrommes  de  repartition  des  erreurs  commises  sur  le 
de  placement  (en  pixel) 


La  methode  du  centre  des  masses  donne  la  repartition  d'erreurla  plus  etroi- 
te,  ce  qui  se  traduit  par  un  faible  ecart  type.  Les  deux  autres  methodes  don- 
nent,  au  contraire,  une  repartition  de  Lerreur  plus  etalee  et  done  un  ecart 
type  plus  important.  Les  erreurs  les  plus  importantes  sont  obtenues  avec  la 
methode  parabolique. 


Tragons  ces  merries  histogrammes  pour  un  rayon  de  1,6  pixel  (fig.  13a  a 
13c).  Le  nombre  total  d'iterations  est  toujours  egal  a  300.  Les  histogram- 
mes  sont,  cette  fois,  beaucoup  plus  pointus  et  done  les  ecarts  types  seront 
beaucoup  plus  petits  que  dans  le  cas  precedent. 
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-0,2  -0,15  -0,1  -0,05  0  0,05  0,1  0,15  0,2 

erreursur  x 
Fehler  aufx 

o)  methode  d'interpolotion :  centre  des  mosses 
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-0,2  -0,15  -0,1  -0,05  0  0,05  0,1  0,15  0,2 

erreursur  x 
Fehler  aufx 


b)  methode  d'interpolotion :  pomboie 


Fig.  13-  Histogmmes  de  repartition  des  erreurs  commises  sur  le 
deplacement  (en  pixel) 


erreur  sur  x 


Fehler  aufx 

c)  methode  d'interpolotion :  goussienne 


4.2.4  Influence  du  deplacement  fractionnaire  sur  I'erreur  de  depla¬ 
cement 

Dans  la  presente  configuration  de  particules,  I'erreur  de  deplacement  ne 
depend  pas  du  deplacement  entier,  mais  uniquement  du  deplacement  frac¬ 
tionnaire  defini  dans  [1  6]  et  rappele  au  paragraphe  3.2.  Pour  etudier  I'in- 
fluence  de  ce  deplacement  fractionnaire,  tragons  pour  differents  rayons  de 
particule  (0,8,  1 ,2  et  1 ,6  pixel)  les  graphes  de  I'erreur  de  deplacement  en 
moyenne  et  en  ecart  type  (fig.  1  4a  a  1  6b)  et  comparons  les  trois  methodes 
d'interpolation.  Les  resultats  ont  ete  obtenus  avec  1  00  simulations  par  va- 
leur  de  deplacement  fractionnaire  (pour  chaque  simulation,  la  position  de 
la  particule  a  Pinstantt]  est  aleatoire). 


40 


ISL-R  119/98 


erreur  moyenne  .  erreur  moyenne 

mittlerer  Fehler  .  ccT  mittlerer  Fehler 


a)  erreur  moyenne  b)  ecort  type 

.14-  Comporoison  des  erreurs  (en  pixel)  obtenues  por  chacune  des  trois  methodes  d'interpolotion  en  fonction  du  deplacement  fractionnaire  pour 
une  particule  de  0,8  pixel  de  rayon 


o)  erreur  moyenne 


b)  ecort  type 


Fig.  15  -  Comporoison  des  erreurs  (en  pixel)  obtenues  por  chacune  des  trois  methodes  d'interpolotion  en  fonction  du  deplacement  fractionnaire  pour 
une  particule  de  1,2  pixel  de  royon 
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emplacement  fractionnaire 
anteilige  Verschiebung 


o)  erreur  moyenne 


0  0,2  0,4  0,6  0,8  1 


emplacement  fractionnaire 
anteilige  Verschiebung 

b)  ecart  type 


Fig.  1 6  -  Comparison  des  erreurs  (en  pixel)  obtenues  par  cbocune  des  trois  methodes  d'interpolation  en  fonction  du  deplacement  fractionnaire  pour 
une  particule  de  1,6  pixel  de  rayon 


Ces  figures  appellent  quelques  commentaires  : 

•  toutes  les  courbes  sont  presque  parfaitement  symetriques  ;  elles  le 
sont  d'autant  plus  que  le  rayon  de  la  particule  est  important ; 

•  quel  que  soil  le  rayon  de  la  particule,  les  courbes  d'erreur  moyenne 
correspondent  aux  trois  methodes  de  determination  sub-pixel  se  re¬ 
couped  pour  un  deplacement  fractionnaire  nul  ou  egal  a  0,5  pixel  : 
la  valeur  prise  par  I'erreur  est  alors  nulle  ; 

•  pour  les  methodes  d'interpolation  parabolique  et  gaussienne,  I'erreur 
moyenne  diminue  lorsque  le  rayon  augmente  ;  ce  sens  de  variation 
est  inverse  pour  la  methode  du  centre  des  masses  ; 

•  les  ecarts  types  sont  d'autant  plus  faibles  que  la  particule  est  de  faille 
importante,  sauf  pour  un  deplacement  fractionnaire  egal  a  0,5  pixel 
pour  la  methode  du  centre  des  masses  ;  dans  ce  cas,  le  pic  de  I'ecart 
type  est  de  plus  en  plus  etroit  et  de  plus  en  plus  haut  ; 

•  ces  figures  sont  en  accord  avec  les  resultats  publies  dans  [1  8]  a  [21  ], 
resultats  obtenus  a  partir  d'autres  approches. 

4.2.5  Influence  du  rayon  de  la  particule  sur  I'erreur  de  deplacement 

Tragons  maintenant  les  courbes  d'erreur  de  deplacement  (en  moyenne  et 
en  ecart  type)  obtenues  avec  les  trois  methodes  d'interpolation  en  fonction 
du  rayon  de  la  particule  lorsque  celui-ci  varie  entre  0,2  et  3  pixel  (fig.  1  7a 
et  1  7b).  Les  variables  du  probleme  sont  la  position  d  I'instant  t]  de  la  parti¬ 
cule  et  son  deplacement.  Afin  de  cadrer  le  couple  de  particules  dans  la  fe- 
netre  d'interrogation  et  de  prendre  en  compte  une  large  gamme  de 
deplacements  dans  la  fenetre  de  32  pixels  x  32  pixels,  les  deux  variables 
seront  generees  comme  suit  : 
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I'abscisse  et  I'ordonnee  de  depart  sont  comprises  entre  1  et  5  pixels, 
les  deplacements  suivant  les  deux  axes  sont  compris  entre  2  et  25 
pixels. 
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Rappelons  que  chaque  point  de  calcul  necessite  1  00  simulations. 


rayon  de  la  particule  (pixel) 
Partikelradius  (Pixel) 

o)  erreur  moyenne 


rayon  de  la  particule  (pixel) 
Partikelradius  (Pixel) 


b)  ecart  type 


Fig.  17  -  Erreurs  de  deplacement  en  fonction  du  rayon  de  la  particule  pour  les  trois  methodes  d'interpolation  sub-pixel 


Notons  tout  d'abord  qu'il  est  impossible  d'utiliser  la  methode  d'interpola¬ 
tion  par  une  gaussienne  pour  des  petits  rayons  (inferieurs  a  0,6  pixel).  En 
effet,  le  pic  de  correlation  de  la  particule  s'etend  alors  sur  moins  de  2  pixels, 
ce  qui  implique  que  Pun  des  niveaux  de  gris  R]  ou  R_]  definis  sur  la  figure 
24  de  [1  6]  est  nul.  Le  calcul  d'interpolation  sub-pixel  etant  obtenu  par  le  lo- 
garithme  des  trois  valeurs  R.-j,  Rq  et  R],  le  resultat  est  complexe  et  ne  cor¬ 
respond  plus  a  une  grandeur  physique.  Ainsi,  la  courbe  correspondent  a 
cette  methode  ne  debutera  que  pour  un  rayon  de  0,6  pixel. 

Notons  aussi  que  les  trois  methodes  d'interpolation  donnent  des  courbes 
semblables  et  qu'elles  possedent  un  minimum  pour  un  rayon  compris  entre 
0,5  et  1  pixel  environ.  II  faudra  done  s'arranger,  lors  de  mesures  de  PIV, 
pour  que  les  particules  qui  sortent  du  generateur  aient  un  rayon  tel  que  le 
rayon  de  leur  image  soit  compris  entre  0,5  et  1  pixel.  Ce  resultat  est,  la  en¬ 
core,  en  accord  avec  les  resultats  publies  dans  [18]  a  [21]. 

Enfin,  signalons  que  des  resultats  concernant  cette  configuration  de  parti¬ 
cules  peuvent  egalement  etre  trouves  dans  [22,  23]. 

4.3  Cas  de  plusieurs  particules  presentes  dans  I'aire 
d'interrogation 

Interessons-nous  maintenant  au  cas  plus  realiste  de  plusieurs  particules  pre¬ 
sentes  dans  la  zone  d'interrogation.  Nous  supposerons  encore  qu'aucun 
bruit  de  fond  ne  perturbe  les  prises  de  vues.  Avant  d'effectuer  les  simula¬ 
tions,  il  faut  preciser  les  conditions  dans  lesquelles  ces  simulations  vont  etre 
realisees. 

4.3.1  Conditions  des  simulations 

Rappelons  qu'en  PIV,  le  nombre  optimal  de  couples  de  particules  dans  une 
zone  d'interrogation  se  situe  autour  de  1  0  (§  4.1 .2).  Nous  essaierons  done 
de  faire  en  sorte  d'avoir  plus  de  5  particules  (pour  que  les  pics  de  correlation 
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soient  suffisamment  eleves)  et  moins  de  15  particules  (pour  eviter  d'even- 
tuelles  superpositions  des  particules).  Contrairement  aux  simulations  du  pa¬ 
ragraphs  4.2,  nous  supposerons  ici  que  des  particules  entrent  dans  la  zone 
d'interrogation  et  en  sortent  entre  les  deux  instants  de  prises  de  vues,  ce  qui 
oblige  a  considerer  la  position  des  particules  dans  une  zone  plus  vaste  que 
la  zone  d'interrogation. 


La  figure  18  represente  un  exemple  de  repartition  de  particules  dans  une 
fenetre  de  64  X  64  pixels  autour  de  la  zone  d'interrogation  de  32  x  32  pixels 
qui  nous  interesse  : 


Fig.  18  -  Position  des  particules  dans  une  fenetre  de  64  x  64  pixels  autour  de 
Poire  dJ interrogation  de  32  x  32  pixels :  la  croix  correspond  a  la 
position  des  particules  a  P instant  tj  et  le  road  a  Pinstant  t2 


Le  nombre  de  couples  de  particules  presents  dans  la  fenetre  de  taille  64  X 
64  pixels  a  ete  initialement  fixe  a  65.  La  repartition  y  etant  aleatoire,  le  nom¬ 
bre  de  couples  de  particules  presents  dans  la  fenetre  de  taille  32  x  32  pixels 
sera  aleatoire  egalement.  La  figure  1  9  represente  le  nombre  de  couples  de 
particules  presents  dans  la  zone  de  32  X  32  pixels  aux  instants  t]  et  \2 /  apres 
1000  simulations  semblables  a  celle  de  la  figure  18.  On  obtient  bien  une 
moyenne  d'une  dizaine  de  couples  de  particules  par  aire  d'interrogation. 


Fig.  1 9  -  Histogramme  de  repartition  du  nombre  de  couples  de  particules  presents 
dans  la  fenetre  d'interrogation  aux  instants  t-j  et  t2 


Zahl  von  Partikelpaaren 


Les  conditions  optimales  de  prises  de  vues  sont  egalement  obtenues  pour 
un  deplacement  total  ^Ax2  +  Ay2  voisin  du  quart  de  la  taille  de  la  zone 
d'interrogation,  done  proche  de  8  a  10  pixels  pour  une  fenetre  de  taille  32 
x  32  pixels  (§  4.1 .2).  Ceci  peut  etre  realise  en  utilisant  pour  Ax  et  Ay,  une 
distribution  normale  de  moyenne  7  pixels  et  d'ecart  type  1  pixel  :  la  figure 
20  montre  I’histogramme  de  repartition  du  deplacement  ^/Ax2  +  Ay2  obte- 
nu  apres  1  000  simulations. 
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Fig.  20  -  Histogromme  de  repartition  du 
deplacement  Jax*  +  Ay2 


4.3.2  Convergence  de  I'erreur  de  deplacement 

Interessons-nous  a  la  convergence  de  I'erreur  de  deplacement  en  moyenne 
et  en  ecart  type.  Pour  cela,  effectuons,  pour  un  rayon  egal  a  1  pixel  et  dans 
les  conditions  du  paragraphe  precedent,  une  serie  de  250  simulations 
(fig.  21a  d  21  f). 


Fig.  21  a, b-  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenue  par  la  methode  du  centre  des  masses  en  autocorrelation  (fig.  21a)  et  en 
correlation  croisee  (fig.  21b) 


d) 


Fig.  21  c,d  ■  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenue  par  la  methode  de  I'ajustement  parabolique  en  autocorrelation  (fig.  21c)  et 
en  correlation  croisee  (fig.  21d) 
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nombre  Alterations  nombre  Alterations 

Zahl  von  Itera  tionen  Zahl  von  Itera  tionen 


Fig .  21  e,f-  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenue  par  la  methode  de  I'ajustement  gaussien  en  autocorrelation  (fig.  21  e)  et 
en  correlation  croisee  (fig.  21  f) 


Les  figures  21  a  a  21  f  correspondant  auxtrois  methodes  d'interpolation  sub¬ 
pixel  considerees  sont  identiques  pour  une  meme  methode  de  correlation. 
On  constate  aussi  que,  contrairement  aux  figures  1  Oa  a  1  1  b,  la  convergen¬ 
ce  n'est  pas  etablie  rapidement,  meme  au  bout  de  250  simulations.  Ceci 
s'explique  simplement  par  les  "sauts11  qui  apparaissent  sur  les  figures  21  a  a 
21  f  et  qui  sont  dus  a  des  erreurs  de  deplacements  importants,  de  I'ordre  de 
plusieurs  pixels  (voir  paragraphe  suivant).  On  remarquera  egalement  que 
les  erreurs  sont  plus  imporiantes  dans  le  cas  de  I'intercorrelation  que  dans 
le  cas  de  ('autocorrelation. 

4.3.3  Distribution  de  I'erreur  de  deplacement 

Les  figures  22a  a  22f  montrent  la  distribution  de  I'erreur  de  deplacement 
dans  le  cas  des  250  simulations  du  paragraphe  precedent.  On  constate 
clairement  que  la  grande  majorite  des  simulations  donne  une  erreur  tres 
proche  de  0  pixel,  alors  que  quelques  simulations  induisent  une  grande  er¬ 
reur  de  plusieurs  pixesl  qui  explique  les  "sauts"  observes  sur  les  figures  21a 
a  2 If,  II  convient  de  preciser  ici  que  sur  les  figures  22a  a  22f,  les  erreurs 
superieures  ou  egales  a  20  pixels  sont  ramenees  sur  le  trait  correspondant 
a  une  erreur  de  20  pixels. 


erreur  commise  erreur  commise 

Fehler  Fehler 


Fig.  22  a,b  ■  Repartition  des  erreurs  de  deplacement  (en  pixel)  obtenues  par  la  methode  du  centre  des  masses  en  autocorrelation  (fig.  22a)  et  en 
correlation  croisee  (fig.  22b) 
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Fig.  22  c,d  -  Repartition  des  erreurs  de  displacement  (en  pixel)  obtenues  par  la  methode  de  I'ajustement  porabolique  en  autocorrelation  (fig.  22c) 
et  en  correlation  croisee  (fig.  22d) 


e) 


erreur  commise 
Fehler 


Fig.  22  e,f*  Repartition  des  erreurs  de  deplacement  (en  pixel)  obtenues  par  la  methode  de  I'ajustement  goussien  en  autocorrelation  (fig.  22e) 
eten  correlation  croisee  (fig.  22f) 


Ces  erreurs  de  plusieurs  pixels  correspondent  a  ce  que  Ton  appelle  couram- 
ment  en  PIV  des  "faux  vecteurs"  et  dont  un  exemple  est  donne  au  paragra- 
phe  suivant.  On  remarquera  alors  que  la  methode  d'autocorrelation  semble 
produire  un  grand  nombre  de  "faux  vecteurs"  donnant  lieu  a  des  erreurs  de 
quelques  pixels,  alors  que  la  methode  de  correlation  croisee  semble  pro¬ 
duire  un  petit  nombre  de  "faux  vecteurs"  donnant  lieu  a  des  erreurs  plus  im- 
portantes  de  I'ordre  de  10  pixels  ou  plus. 

4.3.4  Exemple  de  "faux  vecteur" 

Comme  indique  au  paragraphe  precedent,  il  peut  arriver,  au  cours  d'une 
simulation,  que  le  deplacement  obtenu  differe  de  plusieurs  pixels  du  depla¬ 
cement  impose  aux  particules.  Le  resultat  calcule  correspond  alors  a  un 
"faux  vecteur"  dont  un  exemple  a  ete  obtenu  dans  les  conditions  suivantes 
(la  repartition  de  particules  suivante  est  effectivement  survenue  lors  d'une 
simulation)  : 

•  deplacement  des  particules  :  Ax  =  7,06  pixels  et  Ay  =  8,80  pixels, 

•  a  I'instantt],  24  particules  de  rayon  variable  (moyenne  0,8  pixel). 
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a)  histogramme  suivant  I'axe  x  b)  histogramme  suivont  I'axe  y 


Fig.  V  -  Ajustement ;  par  des  fonctions  parabolique  et  gaussienne >,  des  histogrammes  du  pic  principal  de  la  figure  de  correlation  croisee 
correspondent  a  la  figure  25 


L'ensemble  des  resultats  trouves  lors  du  calcul  des  deplacements  est  repris 
dans  le  tableau  4  ci-dessous  : 


Methode 

Deplacement 
en  x  [pixel] 

Erreur 
en  x  [pixel] 

Deplacement 
en  y  [pixel] 

Erreur 
en  y  [pixel] 

edm,  auto 

10,08 

3,02 

0,06 

-8,74 

edm,  cross 

-1,88 

-8,94 

-19,82 

-28,62 

par,  auto 

10,19 

3,13 

0,12 

-8,68 

par,  cross 

-1,88 

-8,94 

-19,56 

-28,36 

gau,  auto 

10,21 

3,15 

0,15 

-8,65 

gau,  cross 

-1,83 

-8,89 

-19,55 

-28,35 

Tab.  4  -  Resultats  sur  le  calcul  des  deplacements  des  particules 


Rappelons  id  les  abreviations  utilisees  dans  le  tableau  4  : 

•  edm  :  methode  d'interpolation  par  le  centre  des  masses, 

•  par :  methode  d'interpolation  par  ajustement  de  fonction  parabolique, 

•  gau  ;  methode  d'interpolation  par  ajustement  de  fonction  gaussienne, 

•  auto  :  methode  d’autocorrelation, 

•  cross  :  methode  de  correlation  croisee. 

Le  tableau  4  montre  clairement  que  les  deplacements  calcules  sont  tota le¬ 
nient  differents  des  deplacements  imposes.  On  obtient  ainsi  des  erreurs  de 
deplacement  de  plusieurs  pixels  qui  expliquent  les  resultats  des  paragraphes 
4.3.2  et  4.3.3. 

II  apparait  ainsi  egaiement  que  les  faux  vecteurs  ne  sont  pas  uniquement 
dus  a  des  conditions  de  prises  de  vues  des  images  de  PIV  non-optimales 
[1 ,2,1  3],  mais  aussi  a  des  dispositions  aleatoires  particulieres  de  particules 
qui  entrainent  de  fausses  estimations  de  deplacements. 
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4.3.5  Convergence  de  I'erreur  de  deplacement  apres  elimination  des 
"faux  vecteurs" 

Les  "faux  vecteurs"  dont  il  a  ete  question  au  paragraphe  precedent  doivent 
evidemment  etre  elimines  (dans  la  mesure  du  possible)  avant  de  tracer  le 
champ  de  vitesses  resultant.  Differentes  methodes  ont  ete  proposees  et  uti- 
lisees  en  pratique  [1,2,13]  : 

•  comparer  chaque  vecteur  vitesse  avec  ses  voisins  les  plus  proches, 

•  utiliser  un  critere  base  sur  les  hauteurs  respectives  des  pics  de  corre¬ 
lations  les  plus  eleves, 

•  comparer  les  pics  de  correlations  les  plus  eleves  correspondant  d  des 
aires  d'interrogation  de  taille  et/ou  de  positions  differentes. 

En  pratique,  s'il  est  relativement  aise  d'eliminer  les  "faux  vecteurs"  qui  cor¬ 
respondent  a  des  erreurs  de  plus  de  5  pixels,  il  est  par  contre  beaucoup  plus 
difficile  d'eliminer  les  "faux  vecteurs"  qui  correspondent  a  des  erreurs  de  I'or- 
dre  du  pixel.  Dans  la  suite,  nous  allons  supposer  que  ces  methodes  d'elimi- 
nation  des  "faux  vecteurs"  ont  ete  appliquees  a  nos  simulations  et  que  nous 
soyons  ainsi  capables  d'eliminer  les  "faux  vecteurs"  qui  correspondent  a  des 
erreurs  de  plus  de  5  pixels.  De  plus,  nous  supposons  que  nous  pouvons  aus- 
si  eliminer  les  "faux  vecteurs"  qui  correspondent  a  des  erreurs  de  I'ordre  du 
pixel  :  pour  cela,  nous  estimons  d'abord  la  valeur  moyenne  et  Pecart  type 
des  deplacements  simules  et  nous  ne  conservons  que  les  erreurs  de  depla¬ 
cement  qui  se  situent  dans  I'intervalle  centre  sur  la  valeur  moyenne  et  de 
longueur  deux  fois  1'ecart  type.  Cette  operation  est,  dans  la  pratique  des 
mesures  de  PIV,  irrealisable,  mais  elle  donne  une  idee  de  la  convergence 
statistique  des  erreurs  de  deplacement  apres  une  elimination  (ici  ideale)  des 
"faux  vecteurs". 

Les  courbes  de  convergence  de  I'erreur  de  deplacement  ainsi  obtenues, 
avec  toujours  des  particules  de  rayon  1  pixel,  sont  donnees  par  les  figures 
28a  d  28f. 

0.3 
0,25- 
0.2  • 

0.15 
0.1 

0,05  ■ 

0  ■ 

0 

Fig .  28  a,b  -  Convergence  de  I'erreur  de  deplacement  (en  pixel)  obtenues  apres  elimination  des  "faux  vecteurs "  par  la  methode  du  centre  des  masses 
en  autocorrelation  (fig.  28a)  eten  correlation  croisee  (fig.  28b) 


nombre  derations  nombre  d’iterations 

Zahl  von  Iterationen  7a/i/  von  Iterationen 
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nombre  d'iterations 
Zahl  von  Iterationen 


Fig.  28  erf-  Convergence  de  I'erreurde  deplacement  (en  pixel)  obtenues  apres  elimination  des  "faux  vecteurs" par  la  metbode  d'ajustementparfonction 
parabolique  en  autocorrelation  (fig.  28c)  eten  correlation  croisee  (fig.  28d) 
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Fig.  28  e,f-  Convergence  de  I'erreurde  deplacement  (en  pixel)  obtenues  apres  elimination  des " faux  vecteurs"  par  la  metbode  d'ajustementparfonction 
goussienne  en  autocorrelation  (fig.  28e)  et  en  correlation  croisee  (fig.  280 


Ces  nouvelles  courbes  de  convergence  presentent  des  erreurs  de  deplace¬ 
ment  beaucoup  plus  faibles  que  les  courbes  des  figures  21a  a  21  f,  puis- 
qu'elles  ne  depassent  pas  quelques  dixiemes  de  pixel,  De  plus,  on  obtient 
une  convergence  correcte  au  bout  d'une  centaine  de  simulations.  Enfin, 
precisons  que  le  trait  horizontal  sur  ces  figures  correspond  a  Perreur  moyen¬ 
ne  de  deplacement  apres  250  simulations. 

4.3.6  Performances  comparees  des  methodes  de  correlation  et 
d'interpolation  sub-pixel 

Une  autre  fagon  de  comparer  les  performances  des  methodes  de  correla¬ 
tion  et  d'interpolation  sub-pixel  est  de  tracer  la  courbe  representant  le  pour- 
centage  de  simulations  donnant  une  erreur  superieure  a  une  erreur  E0  en 
fonction  de  ce  seuil  E0.  De  ce  point  de  vue,  la  meilleure  methode  est  celle 
qui,  pour  une  abscisse  E0  donnee,  aura  la  plus  petite  ordonnee.  Les  figures 
29a  et  29b  montrent  ces  courbes  dans  le  cas  des  250  simulations  des  pa- 
ragraphes  4.3.2  et  4.3.3,  pour  un  rayon  des  particules  de  1  pixel  (les  "faux 
vecteurs"  sont  pris  en  compte  ici). 

Remarquons  sur  ces  figures  que,  quelle  que  soit  la  methode  choisie,  la 
courbe  decroit  rapidement.  Cependant,  avec  les  echelles  de  ces  figures,  il 
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est  difficile  de  dire  quelle  est  la  meilleure  methode.  Par  contre,  nous  voyons 
clairement  que  la  correlation  croisee  donne  de  meilleurs  resultats  que 
('autocorrelation,  et  que  meme  pour  une  valeur  Eo  de  I'ordre  de  5  pixels,  le 
pourcentage  de  simulations  donnant  une  erreur  superieure  a  Eo  n'est  pas 
nul  :  ce  resultat  est  evidemment  lie  a  la  presence  des  "faux  vecteurs". 


a)  autocorrelation  W  correlation  croisee 

Fig .  29  a,b  -  Pourcentages  de  simulations  dont  I'erreur  est  superieure  a  E0 


Pour  mieux  comparer  les  differences  des  methodes  de  correlation  et  d’inter- 
polation  sub-pixel,  faisons  maintenant  varier  Eo  entre  0  et  0,5  pixel  (fig.  29c 
et  29d). 


c)  autocorrelation 


d)  correlation  croisee 


Fig.  29  c,d  -  Pourcentages  de  simulations  dont !' erreur  est  superieure  a  E0 


En  conclusion,  les  figures  29c  et  29d  montrent  que,  quelle  que  soit  la  me¬ 
thode  de  correlation  utilisee,  c'est  la  methode  d'interpolation  par  ajuste- 
ment  d'une  fonction  gaussienne  qui  donne  les  meilleurs  resultats. 

4.3.7  Influence  du  deplacement  fractionnaire  sur  I'erreur  de  depla¬ 
cement 

Comme  dans  le  cas  d'une  seule  parHcule  presente  dans  I'aire  d'interroga- 
tion  (§  4.2.4),  nous  pouvons  etudier  I'influence  du  deplacement  fractionnai¬ 
re  des  particules  sur  I'erreur  de  deplacement  (fig.  30a  a  30d). 
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Fig.  30a  -  Erm  de  deplacement  en  fonction  da  deplacement 

fractionnoire  (casdel'autocorrelation)  :moyenne  (en  pixel) 
de  I'erreur  de  deplacement 


Fig.  30b  -  Frreur  de  deplacement  en  fonction  du  deplacement 

fraction  naif e  (cas  de  I 'autocorrelation) :  ecart  type  (en  pixel) 
de  I'erreur  de  deplacement 


Fig.  30c  -  Frreur  de  deplacement  en  fonction  du  deplacement 

fractionnoire  (cas  de  la  correlation  croisee) :  moyenne  (en 
pixel)  de  I'erreur  de  deplacement 


Fig.  30d  -  Frreur  de  deplacement  en  fonction  du  deplacement 

fractionnoire  (cas  de  la  correlation  croisee) :  ecart  type  (en 
pixel)  de  I'erreur  de  deplacement 


Les  figures  30a  a  30d  ont  ete  obtenues  pour  un  rayon  des  images  de  par- 
ticules  egal  a  0,8  pixel,  et  chaque  point  de  ces  figure  a  ete  determine  apres 
250  simulations  (les  simulations  correspondant  a  des  "faux  vecteurs"  ont  ete 
eliminees). 

Comme  dans  le  cas  du  paragraphe  4.2.4,  les  courbes  d'erreurs  moyennes 
suivent  une  certaine  symetrie  et  se  recoupent  encore  pour  un  deplacement 
fractionnoire  proche  de  0  ou  0,5  pixel.  La  methode  d'interpolation  par  une 
fonction  gaussienne  donne,  la  aussi,  le  plus  souvent,  I'erreur  la  plus  faible. 

4.3.8  Influence  du  rayon  des  particules  sur  I'erreur  de  deplacement 

Comme  au  paragraphe  4.2.5,  il  est  possible  d'examiner  ['influence  du 
rayon  des  particules  sur  I'erreur  de  deplacement. 
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Fig.  31a  -  Erreur  de  deplacement  (moyenne,  en  pixel)  en  fonction  du 
rayon  des  particules  (cas  de  I'autocorrelation) 


Fig.  31b  -  Frreurde  deplacement  (ecarttype,  en  pixel)  en  fonction  du 
rayon  des  particules  (cas  de  I'autocorrelation) 


Fig.  31c-  Erreur  de  deplocement  ( moyenne ,  en  pixel)  en  fonction  du  Fig.  31  d-  Erreur  de  deplacement  ( ecarttype ,  en  pixel)  en  fonction  du 
rayon  des  particules  (cas  de  la  correlation  croisee)  rayon  des  particules  (cas  de  la  correlation  croisee) 


Dans  ce  cas  precis,  apres  chaque  simulation,  le  resultat  sur  I'erreur  de  de- 
placement  est  teste  pour  savoir  s'il  s'agit  d'un  faux  vecteur  (nous  conside- 
rerons  comme  faux  vecteur  tout  vecteur  dont  I'erreur  de  deplacement  est 
superieure  a  5  pixels)  et  ce  resultat  est  elimine  si  tel  est  le  cas.  Les  figures 
31a  a  31  d  montrent  les  courbes  obtenues  apres  1  00  simulations  par  point 
de  calcul,  pour  des  rayons  de  particules  variant  de  0,2  a  1 ,5  pixel  (pour  les 
autres  conditions  du  calcul,  voir  §  4.2.5).  On  observe,  comme  dans  le  cas 
d'une  seule  particule  par  aire  d'interrogation  (§  4.2.5),  que  I'erreur  de  de¬ 
placement  semble  etre  minimale  pour  un  rayon-image  des  particules  situe 
aux  environs  de  1  pixel. 

Le  fait  que  I'erreur  de  deplacement  semble  etre  minimale  pour  un  rayon- 
image  des  particules  situe  aux  environs  de  1  pixel  peut  se  comprendre  aise- 
ment  en  considerant  les  deux  points  suivants  : 

•  Pour  de  grands  rayons  de  particules,  les  images  de  ces  particules  vont 
immanquablement  se  superposer.  Les  figures  32a  et  32b  en  donnent 
un  exemple  pour  des  particules  de  3  pixels  de  rayon.  Dans  ce  cas  ex¬ 
treme,  on  trouvera  un  "faux  vecteur. 
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Fig .  32a  -  Position  des  particules  de  rayon  egal  a  3  pixels  dans  une 
oire  d' interrogation  de  32  x  32  pixels :  la  croix 
correspond  a  la  position  des  particules  a  I'instant  tj  et 
le  rond  a  Pinstant  t2 


Fig .  32b  -  Image  restituee  par  la  camera  CCD 


•  Lorsque  le  rayon  des  particules  est  petit  (0,4  pixel,  par  exemple),  le 
signal  regu  par  la  camera  CCD  est  tres  faible  puisque  ce  signal  est 
obtenu  par  integration  spatiale  sur  chaque  pixel  du  faible  signal  emis 
par  ces  petites  particules.  Une  eventuelle  superposition  d'images  de 
particules  entraine  alors  une  elevation  locale  forte  du  niveau  de  gris, 
ce  qui  se  traduit  par  un  pic  important  sur  la  repartition  tridimension- 
nelle  des  particules  (fig.  33a  et  33b).  Si  plusieurs  tels  pics  (correspon¬ 
dent  a  des  superpositions)  sont  presents  dans  I'aire  d'interrogation,  la 
correlation  la  plus  forte  est  celle  obtenue  a  partir  de  ces  pics,  ce  qui 
conduit  aussi  a  des  erreurs  importantes  sur  le  deplacement,  voire  a 
des  faux  vecteurs. 


-10  0  10  20  30  40 


axe  x  (pixet) 

Fig.  33a  -  Positions  de  particules  de  rayon  0,4  pixel 


axe  x  (pixel) 


Fig.  33h  -  Diveaux  de  gris  correspondent  a  la  figure  33a  (les  pics 
correspondent  pour  la  plupart  a  des  superpositions  de 
particules) 


4.4  Cas  de  plusieurs  particules  avec  presence  de  bruit 

4.4.1  Aire  d'interrogation  reelle  de  PIV 

Jusqu'a  present,  nous  n'avions  considere  que  des  aires  d'interrogation  ob- 
tenues  apres  un  enregistrement  parfait,  sans  aucune  source  de  bruit.  En 
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realite,  lors  de  prises  de  vues  de  PIV,  les  sources  de  bruit  sont  nombreuses  : 
lumiere  parasite,  bruit  du  CCD  et  des  divers  appareils  electroniques,  imper¬ 
fection  des  optiques,  etc.  [1,2,13].  A  titre  d'exemple,  les  figures  34a  et  b 
montrent  des  representations  bi-  et  tridimensionnelle  des  niveaux  de  gris 
d'une  aire  d'interrogation  obtenue  lors  d'une  prise  de  vue  reelle  de  PIV  ob- 
tenue  a  I'lSL  sur  un  profil  oscillant  en  tunnel  hydrodynamique  [1 4,1 5]  : 


On  observe  bien  les  pics  correspondent  aux  positions  des  particules,  mais 
surtout  on  remarque  la  presence  de  bruit  de  fond  sur  la  figure  34b.  Ce  bruit 
devra  maintenant  etre  pris  en  compte  sur  nos  simulations. 

4.4.2  Simulation  d'aires  d'interrogation  avec  bruit 

Modeliser  separement  chaque  source  de  bruit  intervenant  lors  de  prises  de 
vues  de  PIV  est  une  entreprise  pour  le  moins  delicate.  Par  exemple,  une  ca¬ 
mera  CCD  engendre  a  elle  seule  plusieurs  types  de  bruits  [24,25].  L'influen- 
ce  du  bruit  d'une  camera  CCD  sur  la  precision  de  mesure  en  PIV  a,  d'une 
certaine  fagon,  deja  ete  etudiee  dans  [20],  mais  toutes  les  sources  de  bruit 
n'ont  pas  ete  prises  en  compte. 

Dans  la  suite,  nous  considererons  que  I'ensemble  des  sources  de  bruit  peut 
etre  represente  par  une  variable  aleatoire  gaussienne.  Les  figures  35a  a 
35c  montrent  alors  un  exemple  d'une  aire  d'interrogation  simulee  et  obte¬ 
nue  dans  les  conditions  suivantes  : 

•  50  particules  dans  la  fenetre  de  64  x  64  pixels  avec  des  positions  alea- 
toires  a  I'instantt]  (donnees  par  une  distribution  uniforme), 

•  rayon  des  particules  represente  par  une  distribution  normale  de 
moyenne  1,2  pixel  et  d'ecart  type  0,2  pixel, 

•  deplacement  impose  aux  particules  :  Ax  =  7,05  pixels  et 

Ay  =  1 ,76  pixel, 

•  bruit  gaussien  de  moyenne  20  pixels  et  d'ecart  type  1  0  pixels. 
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par  une  representation  continue,  (.'equivalent  monodimensionnel  de  cette 
represententation  continue  est  schematise  sur  la  figure  37a.  Nous  cher- 
chons  alors  a  calculer  la  fonction  d'autocorrelation  suivante  ; 

+oo 

(37)  cp(t)  =  J  f(x)f(x  +  t)dx  . 

Puisque  nous  avons  f(x)  =  A  pour  0  <  x  <  N,  il  est  facile  de  tracer  f(x+t),  puis 
le  produit  f(x)f(x+t)  et  <p(t)  (fig.  37).  On  obtient  finalement  : 

+eo 

(38)  cp(t)  =  J  f(x)f(x  +  t)dx  =A(N-  It  |)  , 
soit  une  figure  de  correlation  de  forme  triangulaire. 

Le  calcul  qui  vient  d'etre  effectue  est  exactement  I'equivalent  monodimen¬ 
sionnel  et  continu  du  calcul  bidimensionnel  et  discret  qui  a  permis  d'obtenir 
la  figure  36a  (se  rappeler  que  la  figure  36a  a  ete  obtenue  a  partir  de  la 
formule  (31 )  remplacee  ici  par  (37)).  On  comprend  ainsi  que  la  forme  pyra- 
midale  des  figures  36a  et  36b  est  due  directement  a  la  methode  de  correla¬ 
tion  utilisee. 


Une  consequence  directe  de  ce  constat  est  que  les  formules  de  correlation 
simples  (31 )  et  (32)  ne  peuvent  etre  utilisees  dans  le  cas  d'aires  d'interroga- 
tion  avec  bruit  de  fond.  En  effet,  les  figures  36a  et  36b  montrent  que  I'on 
obtiendrait,  a  I'aide  des  methodes  d'interpolation  sub-pixel,  des  deplace¬ 
ments  proches  de  0  pixel  (voir  la  position  du  sommet  de  la  pyramide),  soit 
une  erreur  de  deplacement  Pordre  de  1 00  %. 
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Pratiquement,  dans  le  cas  de  prises  de  vues  reelles  de  PIV,  on  utilise  souvent 
le  coefficient  de  correlation  qui,  dans  le  cas  du  coefficient  de  correlation 
croisee,  s'ecrit  [1  6]  : 


ou  et  Jj/designent  les  moyennes  arithmetiques  de  N^-  et  N7^  (respective- 
ment).  D'autres  auteurs  preconisent  des  formules  de  correlations  differentes 
qui  permettent  des  calculs  numeriques  plus  rapides  :  voir,  par  exemple,  la 
fonction  d'erreur  de  correlation  dans  [26]  et  les  methodes  de  correlations 
binaires  ou  on  "binarise"  au  prealable  les  images  de  particules  (i.e.  on  attri- 
bue  a  chaque  pixel  une  valeur  1  ou  0  suivant  que  ce  pixel  appartient  a  une 
image  de  particule  ou  non)  :  voir  [20,27,28],  par  exemple. 

On  remarque  que  la  difference  principale  entre  ces  methodes  et  les  formu¬ 
les  de  correlation  simples  (31 )  et  (32)  utilisees  jusque  la  dans  ce  Rapport  est 
que  ces  methodes  permettent  d'eliminer  le  bruit  de  fond  avant  de  faire  le 
calcul  de  correlation  proprement  dit. 

4.4.4  Elimination  du  bruit  de  fond 

Dans  ce  paragraphe,  nous  n'allons  pas  utiliser  les  methodes  habituellement 
employees  dans  le  cas  de  prises  de  vues  reelles  de  PIV  (coefficient  de  cor¬ 
relation,  fonction  d'erreur  de  correlation,  binarisation  :  voir  la  fin  du  para¬ 
graphe  precedent),  mais  nous  allons  plutot  mettre  en  oeuvre  une  methode 
simple  d'elimination  du  bruit  avant  de  calculer  les  correlations  a  I'aide  des 
relations  (31 )  et  (32). 

Cette  methode  simple  d'elimination  du  bruit  est  une  methode  de  seuillage  : 
elle  consiste  a  calculer  la  moyenne  1 1  et  I'ecart  type  a  des  niveaux  de  gris 
Nj  j  d'une  aire  d'interrogation,  puis  a  effectuer  I'operation  suivante  : 

Nisj  =  0  si  Njj<|a-hka.  (40) 


Dans  la  suite,  k  sera  appele  "nombre  seuil".  Une  premiere  etude  interessan- 
te  est  d'examiner  I'influence  de  ce  nombre  seuil.  Tragons  alors,  pour  un 
nombre  seuil  variant  entre  0  et  5,  la  figure  donnant  le  pourcentage  de  si¬ 
mulations  dont  I'erreur  de  deplacement  est  1  0  fois  inferieure  apres  elimina¬ 
tion  du  bruit  a  I'erreur  de  deplacement  obtenue  avant  elimination  du  bruit 
(fig.  38).  Chaque  point  de  calcul  est  obtenu  a  partir  de  100  iterations,  les 
conditions  de  simulation  etant  celles  precisees  au  paragraphe  4.4.2. 
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Fig.  38  -  Pourcentage  de  simulations  dont 
Perm  de  deplacement  est  1 0  fois 
inferieure  apres  elimination  du  bruit 
a  Perm  de  deplacement  obtenue 
ovant  elimination  du  bruit 


Notons  que  la  methode  d'interpolation  sub-pixel  n'intervient  pas  dans  cette 
figure.  Les  trois  methodes  donnent,  en  effet,  des  resultats  sensiblement  iden- 
tiques.  La  seule  difference  notable  se  situe  au  niveau  du  choix  de  la  metho¬ 
de  de  correlation  :  la  methode  de  correlation  croisee  donne  de  meilleurs 
resultats  que  la  methode  d'autocorrelation,  et  ce  dans  une  gamme  de  nom- 
bre  seuil  plus  grande. 

La  figure  38  montre  aussi  que  ^elimination  du  bruit  n'est  pas  efficace  pour 
toutes  les  valeurs  du  nombre  seuil.  Si  Ton  fixe  arbitrairement  un  pourcenta- 
ge  d'efficacite  de  67  %,  on  s'apergoit  que,  dans  le  cas  de  ('autocorrelation/ 
le  nombre  seuil  k  doit  etre  compris  entre  0,7  et  4  environ.  Ce  point  peut  se 
comprendre  aisement :  si  le  nombre  seuil  esttrop  petit,  ^elimination  du  bruit 
ne  sera  pas  effectuee  surtous  les  pixels,  ce  qui  fera  apparaTtre  de  nouveaux 
pics  qui  ne  correspondent  pas  a  des  particules  mais  a  du  bruit  (fig.  39a  et 
39b).  Si  ce  nombre  seuil  esttrop  grand,  ('elimination  du  bruit  sera  trop  gros- 
siere  et  seuls  certains  pixels  isoles  resteront  (fig.  40a  et  40b),  ce  qui  detruira 
la  correlation  entre  les  deux  instants  de  mesure.  Dans  les  deux  cas,  les  me¬ 
thodes  de  correlation  (31)  et  (32)  donneront  des  resultats  errones. 


Fig.  39a  -  Exemple  de  niveaux  de  gris  avant  elimination  du  bruit 


Fig.  39b  -  Niveaux  de  gris  apres  elimination  du  bruit  avec  un 
nombre  seuil  egal  a  0,1 
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Fig.  40a  -  Fxemple  de  niveoux  de  gris  ovant  elimination  du  bruit 


Fig.  40b  -  Niveaux  de  gris  apres  elimination  du  bruit  avec  un 
nombre  seuil  egal  a  5 


Tragons  egalement  la  moyenne  de  I'erreur  de  deplacement  avant  et  apres 
elimination  du  bruit  (fig.  41a)  et  b)),  en  fonction  du  nombre  seuil,  a  partir 
des  memes  resultats  que  ceux  qui  ont  servi  a  tracer  la  courbe  38. 


Sch  wellen  wert  Sch  we  I  ten  wert 

a)  courbes  obtenues  en  autocorrelation  b)  courbes  obtenues  en  correlation  croisee 


Fig.  4 1  -  Comparison  des  erreurs  de  deplacement  avant  (courbes  reliant  les  croix)  et  apres  elimination  du  bruit  (courbes  reliant  les  ronds) 


Les  courbes  des  figures  41a)  et  b)  ont  ete  obtenues  en  utilisant  la  methode 
d'interpolation  par  le  centre  des  masses.  Les  deux  autres  methodes  d'inter- 
polation  sub-pixel  donnent  des  resultats  analogues  et  les  courbes  d'erreurs 
sont  sensiblement  les  memes  que  celles  obtenues  ci-dessus.  Ces  courbes 
permettent  de  valider  la  methode  d'elimination  du  bruit  proposee  et  de  con¬ 
firmer  son  efficacite  pour  une  valeur  du  nombre  seuil  comprise  entre  0,8  et 
3,5. 

4.4.5  Influence  du  niveau  de  bruit  de  fond  sur  I'erreur  de  deplace¬ 
ment 

Etudions  maintenant  ['influence  du  niveau  de  bruit  de  fond  sur  I'erreur  de 
deplacement  (fig. 42a  a  43b).  Les  resultats  ont  ete  obtenus  dans  les  condi¬ 
tions  suivantes  : 
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•  deplacement  impose  aux  particules  :  repartition  normale  de  moyenne 
7  pixels  et  d'ecart  type  2  pixels  suivant  les  axes  x  et  y, 

•  rayon  des  particules  :  repartition  normale  de  moyenne  0,8  pixel  et 
d'ecart  type  0,2  pixel, 

•  nombre  seuil  pour  I'elimination  du  bruit  egal  a  2, 

•  200  simulations  pour  chaque  point  des  courbes. 

Les  figures  42a  a  43b  represented,  pour  les  trois  methodes  d'interpolation 
sub-pixel  considerees,  le  pourcentage  des  simulations  donnant  des  erreurs 
superieures  a  0,2  pixel,  en  fonction  du  niveau  moyen  de  bruit  de  fond.  Les 
courbes  des  figures  42a)  et  b)  ont  ete  obtenues  en  fixant  I'ecart  type  a  1  0 
niveaux  de  gris,  alors  que  les  courbes  des  figures  43a)  et  b)  ont  ete  obte¬ 
nues  pour  un  ecart  type  de  50  niveaux  de  gris.  Sur  les  merries  figures  sont 
aussi  representees  les  courbes  obtenues  sans  elimination  du  bruit  de  fond. 


niveau  moyen  bruit  de  fond  niveau  moyen  bruit  de  fond 

mittlerer  Pegel  des  Rauschhintergrunds  mittlerer  Pegel  des  Rauschhintergrunds 


o)  courbes  obtenues  en  autocorrelation  b)  courbes  obtenues  en  correlation  croisee 

Fig.  42  -  Pourcentages  des  simulations  donnant  des  erreurs  superieures  a  0,2  pixel  pour  un  ecart  type  du  bruit  de  fond  egal  a  1 0  niveaux  de  gris 
avant  elimination  (courbes  du  baut)  et  opres  elimination  (courbes  du  bas)  du  bruit 


niveau  moyen  bruit  de  fond  niveau  moyen  bruit  de  fond 

mittlerer  Pegel  des  Rauschhintergrunds  mittlerer  Pegel  des  Rauschhintergrunds 


a)  courbes  obtenues  en  autocorrelation  b)  courbes  obtenues  en  correlation  croisee 

Fig.  43  -  Pourcentages  des  simulations  donnant  des  erreurs  superieures  a  0,2  pixel  pour  un  ecart  type  du  bruit  de  fond  egal  a  50  niveaux  de  gris 
avant  elimination  (courbes  du  haut)  et  opres  eliminab'on  (courbes  du  bas )  du  bruit 
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Ces  figures  appellent  les  commentaires  suivants  : 

a)  la  necessite  de  ('elimination  du  bruit  est  mise  en  evidence.  Sans  cette 
elimination,  il  suffit  d'un  bruit  de  fond  faible  pour  obtenir  pres  de 
100  %  des  simulations  qui  donnent  des  erreurs  toujours  superieures 
a  0,2  pixel, 

b)  le  pourcentage  de  simulations  donnant  des  erreurs  superieures  a  0,2 
pixel  augmente  avec  le  niveau  moyen  de  niveau  de  bruit  de  fond, 
mais  aussi  avec  recart  type  :  en  effet,  a  niveau  de  bruit  moyen  egal, 
le  pourcentage  est  plus  important  lorsque  I'ecart  type  vaut  50  que 
lorsqu'il  vaut  1  0, 

c)  le  choix  de  la  methode  de  correlation  est  a  nouveau  favorable  a  la 
correlation  croisee,  celle-ci  donnant  des  pourcentages  plus  faibles 
par  rapport  a  ^autocorrelation,  a  niveau  de  bruit  egal. 

4.4.6  Conclusion 

Les  differentes  courbes  tracees  dans  ce  paragraphe  montrent  que  le  bruit 
contribue  largement  a  Laugmentation  de  I'erreur  sur  le  deplacement.  II  faut 
done  le  reduire  au  maximum,  d'une  part  en  utilisant  du  materiel  tres  perfor- 
mant  (camera  CCD  generant  un  bruit  de  fond  minimal)  et,  d'autre  part,  en 
s'arrangeant  pour  eliminer  toute  source  lumineuse  parasite.  Toutes  ces  dis¬ 
positions  ne  permettent  cependant  pas  d'eliminer  le  bruit.  II  sera  done  ne- 
cessaire  de  traiter  le  signal  issu  de  la  camera  CCD  en  le  seuillant,  par 
exemple,  suivant  la  methode  proposee  dans  ce  paragraphe. 

4.5  Cas  de  plusieurs  particules  avec  un  gradient  de 
deplacement 

4.5.1  Configuration 

Dans  ce  chapitre,  nous  considerons  le  cas  d'un  ecoulement  avec  un  gra¬ 
dient  de  deplacement.  Par  souci  de  simplicity,  nous  abandonnons  cepen¬ 
dant  le  bruit  de  fond  pour  effectuer  ces  tests.  La  configuration  est  alors  la 
suivante  : 

•  65  particules  presentes  dans  la  fenetre  de  64  X  64  pixels, 

•  rayon  des  images  des  particules  de  1  pixel. 

Pour  simplifier  la  configuration  nous  n'imposerons,  pour  commencer, 
qu'un  gradient  suivant  la  direction  x.  Les  deplacements  imposes  aux  parti¬ 
cules  seront  alors  de  la  forme  : 

AX=  AX0  +  ^^(y- y0),  (40) 

dy 

AY  =  AY0  =  constante  •  (41) 


d(AX)  .  ,  . 

Le  gradient  — ~ —  s  exprime  en  pixel  suivant  x  par  pixel  suivant  y  :  dans  la 

suite,  nous  ne  repreciserons  plus  son  unite.  Le  deplacement  "de  base"  AXg 

correspond  au  deplacement  qu'aurait  une  particule  qui  se  trouverait  a  I'or- 
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donnee  yo  ==16  pixels  lors  de  la  premiere  des  deux  expositions  de  PIV.  C'est 
egalement  ce  deplacement  qui  nous  servira  de  reference  pour  determiner 
I'erreur  commise  sur  le  deplacement.  La  figure  44  montre  une  schematisa- 
tion  des  deplacements  imposes  aux  particules. 


40 


30 


20 


10 

Fig.  44  -  Aire  d' interrogation  de  32  x  32  pixels  montrant  les  deplacements  imposes 

aux  particules  a  cause  du  gradient  d (Ax) /dy.  Le  vecteur  central  &  I'ordonnee  0 
Yq  =16  pixels  represente  AX0 
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4.5.2  Exemple  de  repartition  des  particules 

Tragons  un  exemple  de  repartition  bidimensionnelle  des  particules  pour  la 
configuration  definie  ci-dessus  et  dont  le  gradient  ~  a  pour  valeur 
0,05.  Les  deplacements  “de  base11  suivant  les  deux  directions  ont  pour  va¬ 
leur  AXq  =  8,1 6  pixels  et  y q  =  7,52  pixels  (fig.  45). 


Fig.  45  -  Exemple  de  repartition  des  partkules  correspondent  a  la  configuration 
avec  gradient  de  displacement  d  ( AX) /dy=  0,05 :  la  croix  correspond 
a  la  position  des  particules  a  I'instant  et  le  rond  a  I'instant  f2 
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4.5.3  Influence  du  gradient  de  deplacement  sur  les  pics  de  correla¬ 
tion 

Dans  un  premier  temps,  comparons  les  histogrammes  des  pics  de  correla¬ 
tion  obtenus  avec  gradient  de  deplacement  (voir  equations  40  et  4 1 )  et  sans 
gradient  de  deplacement,  et  ce  pour  la  meme  repartition  de  particules  a 
I'instant  t]  que  sur  la  figure  45  (dans  le  cas  sans  gradient  de  deplacement, 
les  deplacements  imposes  aux  particules  sont  de  la  forme  AX  =  AXq  et 
AY  =  AYq)  :  voir  figures  46a  a  46d. 

Les  deplacements  des  particules  n'etant  pas  tous  les  memes  dans  la  fenetre 
d'interrogation  du  fait  de  ['introduction  du  gradient  de  deplacement,  nous 
nous  attendrions  a  ce  que  le  pic  de  la  fonction  d'autocorrelation  soit  plus 
large  et  moins  haut  que  celui  obtenu  sans  gradient.  Ceci  n'apparaTt  cepen- 
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dant  pas  de  maniere  flagrante  sur  les  figures  46a  a  46d.  Nous  verrons  au 
paragraphe  4.5.6  que  cet  effet  d'elargissement  des  pics  de  correlation  est 
nettement  mieux  mis  en  evidence  pour  des  valeurs  plus  importantes  du  gra¬ 
dient  de  deplacement. 

Remarquons  egalement  que  les  figures  46c  et  46d  sont  tres  semblables,  ce 
qui  etait  attendu  puisque  les  deplacements  suivant  I'axe  y  sont  les  memes 
pour  ces  deux  figures. 


a)  histogramme  obtenu  avec  presence  d'un  gradient  de  deplacement  b)  histogramme  obtenu  sans  gradient  de  deplacement 


Fig.  46  a,b  -  Comparoison  des  histogrammes  suivant  I'axe  x  des  niveaux  de  gris  du  pic  principal  de  la  fonction  d'autocorrelation 


axe  y  axe  y 

c)  histogramme  obtenu  avec  presence  d'un  gradient  de  deplacement  d)  histogramme  obtenu  sans  gradient  de  deplacement 

Fig.  46  c,d  -  Comparaison  des  histogrammes  suivant  I'axe  y  des  niveaux  de  gris  du  pic  principal  de  la  fonction  d'autocorrelation 


4.5.4  Influence  du  gradient  de  deplacement  sur  les  erreurs  de  depla¬ 
cement 

Les  tableaux  5  et  6  montrent  les  deplacements  calcules  par  la  methode 
d'autocorrelation  et  par  les  trois  methodes  d'interpolation  sub-pixel  pour  la 
configuration  de  particules  du  paragraphe  4.5.1  (tableau  5  :  cas  avec  gra¬ 
dient  de  deplacement,  et  tableau  6  :  cas  sans  gradient  de  deplacement). 
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AX 

AXq 

AY 

AY0 

centre  des  masses 

8,09 

8,16 

7,85 

7,52 

parabole 

8,20 

8,16 

7,58 

7,52 

gaussienne 

8,23 

8,16 

7,56 

7,52 

Tab.  5  -  Deplacements  (colonnes  AX  et  A  Y)  cables  par  la  methode  d' autocorrelation  etpar  les  trois  methodes  d'inter- 
polation  sub-pixel  dans  le  cas  d'un  gradient  d(AX)/dy  egal  a  0,05 


AX 

O 

X 

< 

AY 

AYq 

centre  des  masses 

8,18 

8,16 

7,84 

7,52 

parabole 

8,31 

8,16 

7,52 

gaussienne 

8,35 

8,16 

7,60 

7,52 

Tab.  6  -  Deplacements  (colonnes  AX  et  A  Y)  cables  par  la  methode  d'autocorrelation  et  par  les  trois  methodes  dJ inter¬ 
polation  sub-pixel  dans  le  cas  sans  gradient  de  deplacement 


On  remarque,  sur  les  tableaux  5  et  6,  que  les  emplacements  calcules  sont 
beaucoup  plus  proches  des  deplacements  "de  base"  imposes  AXq  et  AYq 
dans  le  cas  sans  gradient  de  deplacement  que  dans  le  cas  avec  gradient  de 
deplacement.  Ainsi,  conformement  a  ce  que  Ton  pouvait  attendre,  un  gra¬ 
dient  de  deplacement  biaise  i'estimation  du  deplacement  des  particules. 

4.5.5  Autre  exemple  de  simulation 

Effectuons  une  deuxieme  simulation  en  prenant,  cette  fois,  un  gradient  plus 
important.  Fixons,  par  exemple,  d(AX)/dy  a  0,1  et  tragons  a  nouveau  un 
exemple  de  repartition  bidimensionnelle  des  particules  (fig.  47)  et  les  histo- 
grammes  des  pics  de  correlation,  mais  uniquement  suivant  I'axe  x  (fig. 48a 
et  b)  puisque  le  gradient  "n'agit"  que  suivant  cette  direction  (les  autres  con¬ 
ditions  de  cette  simulation  sont  celles  des  paragraphes  4.5.1  et  4.5.2). 

Les  tableaux  7  et  8  montrent  ensuite  les  deplacements  calcules  par  la  me¬ 
thode  d'autocorrelation  et  par  les  trois  methodes  d'interpolation  sub-pixel 
pour  la  configuration  de  particules  de  ce  paragraphe  (tableau  7  :  cas  avec 
gradient  de  deplacement,  et  tableau  8  :  cas  sans  gradient  de  deplacement). 


Fig.  47  -  Exemple  de  repartition  des  particules  correspondent  a  la  configuration  avec  gradient 
de  deplacement  d(AX)/dy  =  0,1 :  la  croix  correspond  a  la  position  des  particules 
a  I'instant  tj  et  le  rond  a  /' instant  t2 
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a)  histogramme  obtenu  avec  presence  d'un  gradient  de  deplacement 


b)  histogramme  obtenu  sans  gradient  de  deplacement 


Fig.  48  -  Comparaison  des  histogrammes  suivantl'axe  x  des  niveaux  de  gris  du  pic  principal  de  la  fonction  d'autocorrelation 


AX 

AX0 

AY 

AY0 

centre  des  masses 

7,93 

8,16 

7,16 

7,52 

parabole 

7,80 

8,16 

7,30 

7,52 

gaussienne 

7,78 

8,16 

7,34 

7,52 

Tab.  7  -  Deplacements  (colonnes  AX  etAY)  cables  par  la  methode  d'autocorrelation  etpar  les  trois  methodes  d'inter- 
polation  sub-pixel  dans  le  cas  d'un  gradient  d(AX)/dy  egal  a  0, 1 


AX 

AX0 

AY 

AY0 

centre  des  masses 

7,94 

8,16 

7,19 

7,52 

parabole 

7,92 

8,16 

7,39 

7,52 

gaussienne 

7,89 

8,16 

7,42 

7,52 

Tab.  8  -  Deplacements  (colonnes  AX  et  A  Y)  cables  par  la  methode  d'autocorrelation  et par  les  trois  methodes  d'inter- 
polation  sub-pixel  dans  le  cos  sans  gradient  de  deplacement 

En  comparant  les  figures  48a  et  b  on  observe,  cette  fois,  bien  une  diminu¬ 
tion  de  I'amplitude  et  un  etalement  des  pics  de  correlation  dans  le  cas  d!un 
gradient  de  deplacement.  On  remarque  egalement  que  les  erreurs  sur  I'es- 
timation  du  deplacement  des  particules  ont  augmente  par  rapport  au  cas 
du  gradient  £ga|  ^  0,05  (comparer  les  tableaux  5  et  7). 

dy 

Les  exemples  de  simulations  numeriques  de  ce  paragraphe  ont  permis 
d'avoir  quelques  idees  sur  le  comportement  des  figures  de  correlation  et  de 
I'erreur  de  deplacement  en  presence  de  gradient  de  vitesses.  Des  series  de 
simulations  numeriques  seraient  cependant  encore  necessaires  pour  deter¬ 
miner,  de  fagon  statistique,  ^influence  exacte  des  gradients  de  deplacement 
sur  ('estimation  du  deplacement  des  particules. 
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Dans  ce  Rapport,  les  grandes  lignes  d'un  modele  analytique  de  PIV  video, 
utilise  pour  simuler  la  repartition  des  niveaux  de  gris,  sur  une  matrice  CCD, 
des  images  des  particules  presentes  dans  Pecoulement,  ont  ete  rappelees. 
Ce  modele  analytique  prend  en  compte  les  points  suivants  : 

•  position  ettaille  des  particules, 

•  intensite  du  plan  laser, 

•  fonction  de  transfert  de  modulation  (FTM)  de  Toptique  utilisee, 

•  effets  d'integration,  de  gain  et  de  quantification  de  la  CCD. 

Les  methodes  classiques  de  depouillement  de  telles  images  (autocorrela- 
tion,  correlation  croisee  et  interpolation  sub-pixel)  ont,  elles  aussi,  ete  rap¬ 
pelees.  A  partir  de  ce  modele  analytique,  des  simulations  numeriques  ont 
permis  d'obtenir  des  resultats  statistiques  importants  sur  les  erreurs  commi- 
ses  dans  I'estimation  des  deplacements  des  particules,  et  ce  pour  de  nom- 
breuses  configurations  de  particules  : 

•  une  seule  particule  par  aire  d'interrogation,  sans  bruit  de  fond, 

•  plusieurs  particules  par  aire  d'interrogation,  sans  bruit  de  fond, 

•  plusieurs  particules  par  aire  d'interrogation,  avec  bruit  de  fond, 

•  plusieurs  particules  par  aire  d'interrogation,  avec  gradient  de  depla¬ 

cement. 

Parmi  les  resultats  les  plus  importants  obtenus,  citons  simplement  la  mise  en 
evidence  de  I'apparition  inevitable  de  "faux  vecteurs",  la  determination  du 
diametre  optimal  des  images  des  particules  (environ  1  pixel)  et  I'influence 
importante  du  deplacement  fractionnaire  sur  I'erreur  de  deplacement  des 
particules. 

Les  modeles  analytique  et  numerique  utilises  dans  ce  Rapport  peuvent,  a 
Pavenir,  etre  utilises  pour  etudier  d'autres  methodes  (de  correlation  ou  d'in- 
terpolation  sub-pixel)  qui  permettraient  d'obtenir  de  meilleures  estimations 
des  deplacements  des  particules  en  PIV. 
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